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ABSTRACT 

We investigate the acceleration of shock waves to relativistic velocities in the outer layers of 
exploding stars. By concentrating the energy of the explosion in the outermost ejecta, such trans- 
relativistic blast waves can serve as the progenitors of gamma-ray bursts (GRBs); in particular, 
the "baryon-loading" problem that plagues many models of GRBs is circumvented. Postshock 
acceleration is effective in boosting the kinetic energy in relativistic ejecta. We present physi- 
cally motivated analytic expressions to describe trans-relativistic blast waves in supernovae, and 
we validate these expressions against numerical simulations of test problems. Investigating the 
effect of stellar structure on mass ejection, we find that relativistic ejecta are enhanced in more 
centrally condensed envelopes — e.g., for radiative envelopes, when the luminosity approaches the 
Eddington limit. Convenient formulae arc presented with which to estimate the production of 
relativistic ejecta from a given progenitor. 

We apply our analytic and numerical methods to a model of SN 1998bw, finding significantly 
enhanced relativistic ejecta compared to previous studies. We propose that GRB 980425 is 
associated with SN 1998bw and may have resulted from an approximately spherical explosion 
producing ~ 10~ 6 Mq of mildly relativistic ejecta with mean Lorentz factor T ~ 2, which then 
interacted with a dense circumstellar wind with mass loss rate ~ few x 1CP 4 Mq yr _1 . A highly 
asymmetric explosion is not required. An extreme model of "hypernova" explosions in massive 
stars is able to account for the energetics and relativistic ejecta velocities required by many of the 
observed cosmological GRBs. However, the most energetic bursts require asymmetric expulsion 
of ejecta, perhaps caused by rotationally flattened progenitors. We present simplified models and 
simulations of explosions resulting from accretion-induced collapse of white dwarfs and phase 
transitions of neutron stars. While we find increased energies in relativistic ejecta compared 
to previous studies, these explosions are unlikely to be observed at cosmological distances with 
current detectors, unless extreme explosion energies and asymmetries are invoked. 

Subject headings: gamma rays: bursts — hydrodynamics — shock waves — relativity — supernovae: 
individual (SN 1998bw) — stars: white dwarf, neutron 

1. Introduction 

More than 30 years after their discovery, 
gamma-ray bursts (GRBs) remain one of the most 
outstanding puzzles in astrophysics. A "relativis- 
tic fireball" model has emerged to explain the 
7-ray emission and successfully predict the ex- 
istence of longer wavelength afterglows. How- 



ever, the identity, or identities, of the central en- 
gine driving such fireballs is literally shrouded 
in uncertainty. At least two major problems 
stalk the proposed GRB models. Firstly, extreme 
isotropic energies are implied by the cosmological 
distances of some bursts, determined from red- 
shifted absorption lines in their afterglows. Sec- 
ondly, the "baryon loading - compactness prob- 



lem" results from the large energies, small sizes 
(implied by short timescale variability) and op- 
tically thin spectra of GRBs and requires the 
source to move towards us highly relativistically. 
This in turn mandates miniscule mass-loading, 
- 5.6 x 10- 5 S 7 , 52 (f/100)- 1 e- 1 Mq, where £ 7 , 52 
is the isotropic 7-ray energy in units of 10 52 ergs, 
f = (1 — /3 2 ) -1 / 2 is the mean Lorentz factor of the 
source mass (/3 is the mean velocity relative to the 
speed of light) and e 7 is the conversion efficiency 
of bulk kinetic energy to 7-rays (see, e.g., Piran 
1999, for a review). 

One mechanism that can naturally overcome 
these hurdles is the acceleration of shocks in a pro- 
genitor, which confines and efficiently taps the en- 
ergy from a central explosion resulting from the 
formation of a compact stellar mass object. Al- 
though the shock is initially non- or mildly rela- 
tivistic, it may accelerate to relativistic velocities 
and channel a fraction of the total explosion en- 
ergy into relativistic ejecta. Shock acceleration 1 
occurs in the presence of a steeply declining den- 
sity gradient, and so only a small fraction of the 
ejecta mass reaches relativistic speeds. 

This mechanism has provided the basis for 
many different models of GRBs, ranging from nor- 
mal supernovae (in the first astrophysical GRB 
model; Colgate 1974), to extremely energetic su- 
pernovae or "hypernovae" (Paczynski 1998), to 
phase transitions of neutron stars into strange 
stars (Fryer & Woosley 1998). Accretion-induced 
collapse of white dwarfs to neutron stars (Woosley 
& Baron 1992; Fryer et al. 1999) may also drive 
a shock into the outer atmosphere of the white 
dwarf, but such shocks have so far been exam- 
ined only in the context of much weaker Type 
la supernova explosions (Berezinsky et al. 1996). 
Shock acceleration is also an important factor in 
aspherical models such as those involving rotation- 
ally flattened progenitors (Chevalier & Soker 1989; 
Khokhlov et al. 1999) and jet-like explosions of the 
core (MacFadyen & Woosley 1999; Hoflich et al. 
1999). Even if a jet pierces the stellar envelope 
entirely (as in the model of Aloy et al. 2000), the 
production of relativistic ejecta by shock acceler- 
ation may dominate the observed GRB in most 

We shall use the term "shock acceleration" to refer to this 
hydrodynamical effect, which ought not be confused with 
particle acceleration in shocks. 



directions. 

While the propagation of shock waves is well 
understood in the nonrelativistic (Sedov 1946, 
1959; Taylor 1950; Ostriker & McKee 1988; 
Matzner & McKee 1999) and ultrarelativistic 
(Johnson & McKee 1971; Blandford & McKee 
1976) regimes, similarity solutions do not exist for 
the mildly relativistic case. Thus, previous ana- 
lytic models of the acceleration of shocks at these 
intermediate velocities have extrapolated non- 
relativistic results with simple scaling laws (e.g. 
Gnatyk 1985) intermediate between the non- and 
ultrarelativistic limits. Furthermore, the post- 
shock evolution of trans-relativistic ejecta has re- 
ceived only a crude treatment, if any, in previous 
work. 

The aim of this paper is to provide a quanti- 
tative, theoretical description of trans-relativistic 
shock and postshock acceleration, which will be 
of use in many GRB models. A trans-relativistic 
blast wave is one in which some of the mass is 
accelerated to highly relativistic velocities, but 
most of the mass is not; in particular the ratio 
of the explosion energy to the rest energy of the 
ejecta is < 0(1). In §2, utilizing a relativistic 
one-dimensional Lagrangian hydrodynamic code 
based on the method of exact Riemann solvers 
(Marti & Miiller 1994), we investigate the dynam- 
ics of shocks and shocked ejecta in idealized model 
density distributions. We present physically mo- 
tivated analytic approximations to describe the 
shock and postshock acceleration as well as the 
kinetic energy distribution of ejecta from density 
distributions approximating polytropes at their 
edges. We explore how a progenitor's structure 
affects its ejecta in §2.4, finding that highly con- 
densed envelopes produce fast or relativistic ejecta 
most efficiently and that this corresponds to stars 
whose luminosities approach the Eddington limit. 
We present formulae to estimate the energy in 
relativistic ejecta from the parameters of a given 
progenitor star. Our results, although primarily 
derived for spherically symmetric explosions, can 
also be applied to sectors of aspherical explosions, 
a topic we address in §2.5. Explosions where the 
effect of gravity is strong, for example relevant to 
ejecta from neutron stars, are treated in §2.6. 

We apply our results to four astrophysical ex- 
amples, as potential GRB engines (§3). First, and 
in most detail, we consider supernova explosions in 
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the exposed cores of massive post main-sequence 
stars (§3.1). The probable association of super- 
nova (SN) 1998bw with GRB 980425 (Galama et 
al. 1998) has sparked renewed interest in such 
models. Other supernovae for which a suspicion 
of GRB association exists include SN 1997cy with 
GRB 970514 (Germany ct al. 2000; Turatto ct al. 
2000) and SN 1999E with GRB 980901 (Thorsett 
& Hogg 1999). The late time excess flux, rela- 
tive to that expected from a power-law decline, 
of the afterglows of GRB 970228 (Fruchter et al. 
1999) and GRB 980326 (Bloom et al. 1999) has 
been cited as evidence for the presence of super- 
novae (Bloom et al. 1999; Reichart 1999; Galama 
et al. 2000). However, alternative mechanisms of 
afterglow re-radiation (Waxman & Drainc 2000) 
and scattering (Esin & Blandford 2000) by dust, 
have also been proposed to produce the excess 
flux. Cosmological GRB afterglows are associated 
with the stellar population of their host galaxies 
(Bloom et al. 2000b), but there is only tentative 
evidence for a direct association with star forming 
regions (Kulkarni et al. 2000; Djorgovski, Bloom, 
& Kulkarni 2000; Galama & Wijers 2000), which 
would be expected if they arose from the deaths 
of massive stars. 

If there is an association between GRB 980425 
and SN 1998bw, then the 7-ray energy is only 
~ 10 48 ergs, thousands to millions of times weaker 
than the inferred isotropic energies of the cosmo- 
logical bursts; moreover, its associated galaxy is 
much too near to be drawn from the same distri- 
bution as the other known host galaxies (Hogg & 
Fruchter 1999). This event may be an example 
of a subclass of GRBs, termed S-GRBs (Bloom 
et al. 1998a). Woosley et al. (1999) simulated an 
energetic (~ few x 10 52 ergs) explosion in a va- 
riety of carbon-oxygen and helium star progen- 
itors, attempting to fit the observed supernova 
lightcurvc. They considered the possibilities of 7- 
ray emission from shock break-out and the interac- 
tion of ejecta with a circumstellar wind. With non- 
relativistic codes they found the simple spherical 
case unable to produce enough relativistic ejecta 
to account for even this weak GRB. This result, 
together with the extremely large isotropic ener- 
gies (Ej ~ 10 51 — 10 54 ergs) of the cosmologi- 
cal bursts, has focussed attention on aspherical 
explosions involving jets, which may be able to 
explain the different types of GRBs as a func- 



tion of the viewing angle to the jet (Wheeler et 
al. 2000a,b; MacFadyen & Woosley 1999; Mac- 
Fadyen et al. 2000; Aloy et al. 2000; Wang et al. 
2000). However, uncertainties in previous models 
of trans-relativistic dynamics warrant further in- 
vestigation of the simplest spherical case. Indeed, 
our relativistic analysis of the model considered 
by Woosley et al. (1999) yields sufficient mildly 
relativistic ejecta to account for GRB 980425 and 
radio observations of SN 1998bw, for a progeni- 
tor embedded in a dense wind with mass loss rate 
<~ few x 10 _4 M Q yr _1 . Our analysis therefore val- 
idates the suggestion by Matzner & McKee (1999) 
that GRB 980425 could be the product of a spher- 
ical explosion of SN 1998bw. 

In §3.2 we continue with explosions in the cores 
of massive stars, but now at much greater ener- 
gies. These hypernovae 2 are hypothesized to result 
from magnetic extraction of the rotational energy 
of nascent neutron stars or black holes (Paczyn- 
ski 1998). We consider an extreme example with 
5 x 10 54 ergs, which might arise from the forma- 
tion of a ~ 10 M Q black hole. For the purposes 
of this crude estimate, we embed the explosion in 
the ~ 5M envelope ejected from our SN 1998bw 
progenitor, and calculate the resulting energy dis- 
tribution. To compare to the observed cosmolog- 
ical bursts we derive conservative constraints on 
the minimum mean Lorentz factor of ejecta carry- 
ing the 7-ray energy, by considering the timescale 
of the burst, and requiring the implied circumstel- 
lar material to be optically thin to the burst pho- 
tons. We also consider the minimum Lorentz fac- 
tor required for the optical depth due to photon- 
photon pair production to be small. We find that 
the overall energetics and relativistic ejecta veloc- 
ities required by many cosmological bursts can be 
explained by such a simple, though admittedly ex- 
treme, model. However there are some events, no- 
tably GRB 990123, that lie well beyond the reach 
of the spherical model. Considering aspherical ex- 
plosions (§2.5), we calculate that enhancements in 
shock strength along the preferred axis, perhaps 
set by rotation, by factors of only a few are re- 
quired to explain all cosmological GRBs. Interac- 
tion of the relativistic ejecta with inhomogeneitics 
may give bursts some short timescale substructure 

2 Unlike Iwamoto et al. (1998), we shall reserve the term "hy- 
pernovae" for much more energetic events than SN 1998bw. 
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(Dermer & Mitman 1999; Fcnimorc et al. 1999), 
though this can only account for a small fraction of 
the total burst energy (Sari & Piran 1997). Insta- 
bilities associated with an accelerating shock (Luo 
& Chevalier 1994) may create postshock velocity 
perturbations that also lead to GRB variability. 

Finally wc consider simple models of explosions 
in white dwarfs and neutron stars (§3.3), in order 
to investigate shock and postshock acceleration in 
their outer layers. Such shocks may be driven by 
the accretion-induced collapse of a white dwarf to 
a neutron star or the phase transition of a neutron 
star to a hypothesized "strange" star, composed of 
more stable strange quark matter. For the produc- 
tion of relativistic ejecta, these models have the 
advantage that large explosion energies may be 
coupled to very small ejecta masses. Given a total 
ejecta mass and energy, as predicted from detailed 
core-collapse models (Fryer & Woosley 1998; Fryer 
et al. 1999), we estimate the energy distribution 
of the ejecta. Again, we find enhanced amounts 
of energy in relativistic ejecta compared to previ- 
ous estimates. We compare to cosmological GRBs, 
this time by combining the burst timescale with an 
assumed density of a uniform ambient medium. 
We find that white dwarf explosions are unlikely 
to be relevant to these bursts. Our neutron star 
models, while more energetic, need to be pushed 
to extreme limits to provide observable fluxes at 
cosmological distances. 

2. Theory 

The process of mass ejection in a supernova 
explosion proceeds in a sequence of distinct dy- 
namical phases (Matzner & McKee 1999). In the 
first phase, the stellar envelope is engulfed by the 
outward-propagating blast wave launched by col- 
lapse of the core. Except perhaps in the vicinity of 
the central engine and within the shock itself, the 
blast wave is adiabatic. For an explosion in the 
tenuous outer envelope of a supernova, the post- 
shock pressure is dominated by photons; for a rela- 
tivistic shock in a compact object, it may be dom- 
inated by relativistic motions of leptons (Weaver 
1976). In either case, the ratio of specific heats 
is 7 = 4/3. The boundary of this blast wave is 
a shock front, which imprints its velocity (and an 
associated specific entropy) on the postshock gas. 
The shock velocity decreases in most parts of the 



star, but it increases in regions with a sufficiently 
steep density profile - an effect analogous to the 
cracking of a whip. As the postshock pressure is 
dominated by radiation pressure, the shock front 
has a finite optical depth (for which photons dif- 
fuse upstream as fast as gas flows downstream). 
Thus the shock breaks out of the star at some fi- 
nite depth. Only for a sufficiently energetic explo- 
sion in a sufficiently compact star will the shock 
accelerate to relativistic velocities before it breaks 
out. In §2.1 we investigate shock acceleration in 
such progenitors. 

After the stellar envelope has been shocked, its 
thermal and kinetic energies are approximately 
equal. The next phase of evolution is one of post- 
shock acceleration, in which heat is converted into 
outward motion and the ejecta approach a state 
of free expansion. For the gas in the outer lay- 
ers of the star, this increases the velocity by over 
a factor of two compared to the postshock state 
if the shock is nonrelativistic (Matzner & McKcc 
1999). However, gas hit by a relativistic shock 
undergoes even more dramatic acceleration, be- 
cause postshock kinetic and thermal energy densi- 
ties are both dominated by photons, which even- 
tually transfer their energy to the baryons in the 
freely expanding state. This conversion of internal 
energy to kinetic leads to a large increase in the 
Lorentz factor during the acceleration phase. The 
acceleration is accentuated by the fact that each 
gas element is pushed outward by interior ejecta, 
which though slower, are at higher pressure (§2.2). 

Each stage of this process has been the study 
of previous theoretical investigations; our goal is 
to investigate in greater detail the connection be- 
tween the shock velocity and the final expansion 
velocity, across the transition from Newtonian to 
relativistic motion. To describe this process ana- 
lytically, we present formulae that have been cali- 
brated against numerical simulations for a variety 
of progenitor stars (§§2.1 and 2.2). In §2.3 we 
give an analytical prediction of the distribution 
of energy with Lorentz factor in the state of free 
expansion, i.e., the information necessary to calcu- 
late the circumstcllar interaction and the resulting 
burst of energetic radiation. We consider the de- 
pendence of relativistic mass ejection on the pro- 
genitor's structure in §2.4. Wc apply our dynam- 
ical results to sectors of aspherical explosions in 
§2.5 and explosions where the ejecta are strongly 
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affected by gravity in §2.6. 



2.1. Shock Propagation 



The propagation of a nonrelativistic blast wave 
through the stellar interior and the subsequent ac- 
celeration phase were studied in detail by Matzner 
& McKee (1999). The shock velocity, (3 s c, re- 
sponds to two competing trends. First is the over- 
all deceleration as mass is entrained through the 
shock, P s oc [_Ei n /m(r)] 1//2 , where 



m{r) 



m(r) - M„ 
Ml; 



0<m<l (1) 



is the fraction of the total ejecta mass, M e j, within 
r (excluding any remnant mass M rom ), and 



E, 



M ci c 2 



(2) 



is the injected energy in units of the rest en- 
ergy, which is < 0(1) for trans-relativistic explo- 
sions. Second is the acceleration of the shock front 
down the declining density gradient, for which 
(3 S oc [rh{r)M e -J{pr z )] a ™ where a nr ~ 0.19 for 
typical outer density distributions (Sakurai 1960). 
The shock propagation is a combination of these 
two trends, and it is not surprising that a product 
of the two formulae yields an excellent approxi- 
mation to the shock velocity (Matzner & McKee 
1999): 



Ps=A 



m 



1/2 



f mM c] 
\ pr 3 



(3) 



The coefficient A can be evaluated for a spheri- 
cal Sedov blast wave in the distribution p oc r~ fep , 
in which it takes the value 



A 



-1/2 



4-7T 

i-k t 



(4) 



where a = E- m / (rhfll), constant for self-similar 
blast waves, measures the fraction of the explosion 
energy that is kinetic. Note that the shock accel- 
erates in this region only if k p > 3. Matzner & 
McKee (1999) found that the Primakoff distribu- 
tion k p = 17/7 (for which A = 0.794) adequately 
describes the interiors of red and blue supergiant 
progenitors of type II supernovae. However, we 
have found that the accuracy of equation (3) is en- 
hanced if A is evaluated for the typical value of k p 
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Fig. 1. — Spherical shock propagation down a den- 
sity profile p — phx 3 , where x = 1 — r/R, for 
different explosion energies, E ln . Simulation re- 
sults {solid lines) are compared to the analytic pre- 
diction of this paper (eqs. [14] and [15]) {dashed 
lines), and to the variation of the scaling predicted 
by Gnatyk (1985) (GMM scaling) as expressed 
by equation (9) {dotted lines). The normalization 
constant, A, is 0.68 for E in < 0.03, and 0.72, 0.74 
for E[ n = 0.1,0.3, respectively. 
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Fig. 2. — Spherical shock propagation down a real- 
istic stellar density profile, given by equation (11). 
The differences between the GMM scaling (eq. [9]) 
(dotted lines) and our improved formula (eqs. [14] 
and [15]) (dashed lines) are clearly revealed in the 
outer region (x 1 / 3 < 0.8), where the density gra- 
dient is very steep. The simulation results are 
again shown by the solid lines. The normaliza- 
tion constant, A, is 0.736 for all energies. The 
lower panel shows the ratio of the analytic pre- 
dictions to the numerical results. For both cases 
E in = 0.0033, 0.03, 0.3 from top to bottom. The 
GMM scaling under predicts relativistic shock ac- 
celeration and over-predicts nonrelativistic. The 
final Lorentz factor of ejecta depend sensitively 
on T s (§2.2), and the observed difference, ~ 33%, 
between the two laws results in a difference of a 
factor of ~ 3 in final Lorentz factors. 



in the interior of progenitor. To make this evalua- 
tion, we use the pressure gradient approximation 
of Ostriker & McKee (1988), with 7 = 4/3: 



2(29-7fc p )(3-fc p ) 
113- 35fc„ : 



(5) 



which is exact for the Primakoff blast wave, k p — 
17/7. For the progenitor model C06 of Woosley 
ct al. (1999), which we use to simulate SN 1998bw 
(§3.1), we estimate k p ~ 1.9 which gives A ~ 
0.736. 

We also consider the velocity normalization 

for shocks that develop in centrally concentrated, 

non-power-law distributions, with monotonically 

steepening density gradients as r — ► R. These 

types of distribution mimic those of white dwarfs 

and neutron stars. We estimate a mean k p in the 

central region, where the shock structure forms, 

by evaluating a moment of the whole density dis- 
t-R 



tribution, J z = j"(r/R) l d(m/M). 
law distribution 



kp — 



(3 + l)Ji - 3 
Ji-l ' 



For a power 



(6) 



From numerical simulation (below), we find I = 
— 1 gives an estimate of an effective k p which ac- 
curately predicts A, via equations (4) and (5), for 
Ei n <C 1 , in several different distributions we have 
tested. We find A may have a weak energy depen- 
dence for Ei n ~ (9(1) (below). 

Although equation (3) is quite accurate in the 
nonrelativistic regime, it breaks down when (3 S — ► 
1 and the shock becomes trans-relativistic. Even 
in the context of a planar, power-law density dis- 
tribution, the dynamics are not self-similar in 
this stage. Therefore, it has received little atten- 
tion in previous work, whereas the opposite limit 
of self-similar, ultrarelativistic shock propagation 
was treated approximately by Johnson & McKee 
(1971). However, a knowledge of the nonrela- 
tivistic and ultrarelativistic shock velocity scal- 
ings does not suffice to specify the relative co- 
efficient connecting the two stages. Moreover, it 
does not allow one to describe ejecta that are only 
mildly relativistic; as we shall show, these trans- 
relativistic ejecta contain more energy than the 
ultrarelativistic. For these reasons, we wish to in- 
vestigate in detail the transition between Newto- 
nian and ultrarelativistic motion. 
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Using an approximate treatment, Johnson & 
McKee (1971) presented the following estimate of 
the shock Lorentz factor, T s = (1 — /3^) -1 / 2 , valid 
in the ultrarelativistic (T s ^> 1) limit: 



r ~ 



( 



Ei 



V2 VPC 5 



(7) 



where a r = ±V3(2 + V^T 1 ~ 0.232, and Ei is the 
energy density at some reference point in the flow. 
Johnson & McKee obtained this result by assum- 
ing that initially a region of hot gas, with a spa- 
tially uniform value of E i} abutted cold gas with 
decreasing density p(r). As a result, the Riemann 
invariant was the same on all forward character- 
istics emanating from the hot gas. In the case 
considered here there is no region with constant 
Ei, so the exact values of a r and the numerical 
coefficient may be somewhat different. 

Noting the similarity between the Newtonian 
and relativistic indices a r and a nr that relate 
the shock velocity to the density profile, Gnatyk 
(1985) suggested the following formula for shocks 
bounding spherical blast waves: 



T s p s oc (pr 3 



(8) 



where a g ~ 0.5 for decelerating and 0.2 for ac- 
celerating shocks. This approximation was found 
to be inferior to equation (3) in the nonrelativistic 
limit by Matzner & McKee (1999), but it does sug- 
gest that the quantity T s f3 s can be approximated 
analytically. 

However, equation (8) does not specify the nor- 
malization of the shock velocity. We therefore con- 
sider instead the following formula: 



1/2 



(r s /3 s )GMM — A 



M ( 



171 



rhM c 



\ pr 3 



(9) 



In the Newtonian limit this formula reduces to a 
shock velocity approximation very similar to equa- 
tion (3), and it makes a transition to relativistic 
motion in the manner suggested by Gnatyk (1985), 
(equation [8]): hence the subscript GMM. Since 
to ~ 1 in the region of strong acceleration at the 
outside of a star, this formula has the same scal- 
ing as equation (8) there. Whereas Woosley et 
al. (1999) chose the shock-velocity coefficient nu- 
merically, equation (9) uses our earlier results and 
those of Matzner & McKee (1999) to prescribe it 



analytically. Note that it is Gnatyk's exponent 
a g = 0.2, rather than a nr = 0.19, that enters this 
expression. Insofar as equation (3) accurately pre- 
dicts the behavior of a nonrelativistic shock front 
with a DI — 0.2 (and we find that it does), equa- 
tion (9) represents the result of an extrapolation to 
relativistic velocities using Gnatyk's scaling. It re- 
mains to be determined whether the transition to 
relativistic motion is captured accurately by this 
formula. 

To study the dynamics of trans-relativistic 
blast waves, we have developed a relativistic, 
one-dimensional Lagrangian hydrodynamics code 
based on an exact Riemann solver (Marti & Miillcr 
1994). The Lagrangian nature of our code is a key 
requirement for handling large ranges in density 
with fine zoning, of particular relevance for simu- 
lating shocks accelerating in stellar atmospheres. 
The code accurately handles standard test prob- 
lems, including Centrella's (1986) proposed wall 
shock and Riemann shock tube tests. Additionally 
the code recovers the ultrarelativistic Blandford- 
McKee (1976) solution for spherical blast waves. 

With this code we have conducted a suite of 
numerical calculations. For much of this study we 
employ a "stellar" model density distribution: an 
outer stellar envelope of polytropic index n has the 
density profile 



p = p h {R/r-l) n , 



(10) 



in any region where the external mass to ox can 
be ignored compared to M* and where n does 
not vary. If this region extends inward to half 
the radius (e.g., for polytropes with n > 3; Chan- 
drasekhar 1939) then the coefficient ph is equal to 
the density there: ph = p(r = R/2). For our stel- 
lar model, we join this outer profile to an inner, 
power-law density distribution: 



Ph x 



1 



r <r c 
r c <r < R. 

(11) 

We choose a core radius r c ~ -R/2, polytropic in- 
dex n = 3.85, and power law index k p = 1.9 to 
resemble the progenitor model for SN 1998bw con- 
sidered in § 3.1. In order to find analytical formu- 
lae valid for various density distributions, we have 
additionally considered an "external power law en- 
velope" model, 

P = p h x n , (12) 



7 



also with n — 3, where 

x = l-^- (13) 

Equation (12) is the x — > limit of the polytropic 
profile (10), but significantly underestimates p in 
the interior (e.g., by a factor of 2™ at x = 1/2). 

After an injection of thermal energy in the cen- 
tral few zones, a shock develops and propagates 
outwards, decelerating in the inner region, where 
the density gradient is shallow, and accelerating 
in the outer layers, where the density gradient is 
steep. We have investigated explosions with E- ln 
ranging from 0.001 to 0.33. In each case, the res- 
olution was increased either until the numerical 
results converged, or until a correction for incom- 
plete convergence could be ascertained. 

Empirically, we find that the scaling proposed 
by Gnatyk (1985), as realized in equation (9), is 
roughly correct but not sufficiently accurate for a 
study of the energetics of relativistic ejecta. In 
order to capture the change in the value of the 
exponent from a nr to a r , we adopt instead the 
following generalization of equation (9): 

r s (3 s = P (i + p 2 r 12 , (w) 

> = -(C(^r <i5> 

For the density distributions described above, our 
numerical results give a nr = 0.187, in agreement 
with the nonrelativistic self-similar theory and 
equation (3). We adopt this numerical value in 
our subsequent formulae. enters in many of 
our expressions and we set this equal to 5.35. In 
the ultrarelativistic regime (p> 1), equation (14) 
was chosen to reproduce the approximate result 
a r ~ 0.232 obtained by Johnson & McKee (1971). 

In Figures 1 and 2 we compare equations (9) 
and (15) to the results of our model simulations of 
explosions in the "external power law" and "stel- 
lar" density distributions. The two shock velocity 
formulae under consideration are extremely sim- 
ilar, differing by only very small powers of their 
parameters in both the non- and ultrarelativistic 
limits. Their differences are most clearly seen in 
the stellar model, where density, and hence shock 
velocity, spans a greater range. 

Limited numerical resolution does not allow us 
to follow the acceleration of a nonrelativistic shock 



into the ultrarelativistic regime. Instead we have 
increased the explosion energy to investigate ac- 
celeration in different velocity regimes. Our re- 
sults show that our improved shock propagation 
model (eq. [14]) is accurate to within ±5% for all 
energies and density distributions tested. 3 Fur- 
thermore we find the simple Gnatyk scaling un- 
der predicts the acceleration of relativistic shocks 
and over-predicts the acceleration of nonrelativis- 
tic shocks. Our highest energy explosion of the 
stellar model shows that the Gnatyk scaling un- 
derestimates, by ~ 33%, the value of T s reached 
by a shock that accelerates from T S (3 S ~ 1 to 
T S (3 S ~ 10. The final Lorentz factor, Tf, achieved 
by shocked ejecta depends sensitively on T s (see 
§2.2). For example, postshock acceleration am- 
plifies the above decrement into a factor of three 
error in Tf. 

We must note that our approximation (eqs. [14] 
& [15]) gives T s oc (E in /m) - 56 rather than the 
correct scaling T s oc (Ein/m) 1 / 2 that obtains in 
spherical, ultrarelativistic blast waves (Blandford 
& McKee 1976). Because of this, the accuracy 
of our formula breaks down for blast waves that 
start out ultrarelativistically; i.e. with E- ln > 1. 
We begin to see these deviations by the time 
E- m ~ 0.3, as the value of A tends to increase 
slightly (< 10%), for a given density distribution. 
Equations (14) and (15) are strictly valid only for 
trans-relativistic blast waves, i.e. for shocks that 
are non- or semi-relativistic to begin with, but that 
may accelerate to high T s , i.e. E m < 0(1). This is 
well satisfied for most of the astrophysical contexts 
to which we shall apply the formula in §3, e.g., 
our model for supernova 1998bw has E m ~ 0.003. 
However, the extreme hypernova models consid- 
ered in §3.2 are at the limits of this approxima- 
tion's validity, with E- m ~ 0.6. 

2.2. Postshock Acceleration 

After the entire stellar envelope has been tra- 
versed by the shock wave, it is left in a state of 
outward motion. However, its heat content is com- 
parable to its kinetic energy, so it will accelerate 
considerably as heat is converted into motion. As- 

This is after allowing for the normalization constant, A, 
to increase slightly with E in for energetic explosions, e.g., 
for the external power law envelope model, the low energy 
(i?i n < 0.03) explosions have A = 0.68 as predicted from 
cq. (6), while for E in = 0.1,0.3, A = 0.72,0.74. 
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suming a strict conservation of lab- frame energy in 
each mass element - from the postshock state to a 
state of free expansion - allows one to roughly esti- 
mate the magnitude of this acceleration. However, 
this is a significant underestimate of the accelera- 
tion since the outermost layers gain energy at the 
expense of the slower ejecta within them, because 
of a steeply declining pressure gradient. 

The relativistic strong shock jump conditions 
give (Blandford & McKee 1976) 

F 2_ (r 2 + i)[7(r 2 -i) + i] 2 

is 7(2-7)(r 2 -l) + 2 ' 1 ' 

where the subscript 2 refers to the fluid immedi- 
ately behind the shock. For 7 = 4/3, we have 
P s /P 2 = 7/6 and T s /r 2 = V2 in the non- and 
ultrarelativistic limits, respectively. 

Now consider the postshock acceleration of ma- 
terial that was struck by a nonrelativistic, strong 
shock as it approached the periphery of the pro- 
genitor star. Thermal and kinetic energies are in 
perfect equipartition behind the shock. There- 
fore, strict energy conservation in each element im- 
plies P f //3 2 = \/2 and so (3 f /0 a = 6^2/7 = 1.21. 
This should be contrasted with the results of self- 
similar shock acceleration and postshock acceler- 
ation, which for 7 = 4/3 yield Pf//3 S = 2.16 and 
2.04 for n = 3/2 and 3, respectively (Matzner & 
McKee 1999). Thus, the energy per baryon in- 
creases by factors of 3.2 and 2.8 in these two cases 
due to the concentration of energy in the fastest 
ejecta. 

This phenomenon is even more dramatic in 
the context of relativistic acceleration, where the 
postshock energy density is dominated by inter- 
nal energy rather than the rest mass of the fluid. 
For a strong relativistic shock, the total post- 
shock energy per unit rest mass, evaluated in the 
fixed frame and neglecting terms of order r,T 2 , 
is 41^/3 = 2r 2 /3 (Blandford & McKee 1976). 
If each baryon maintained exactly this share of 
energy into the final state of cold, free expan- 
sion, then 2r 2 /3 would also equal its terminal 
Lorentz factor Tf. The theory of acceleration 
behind an ultrarelativistic shock was treated ap- 
proximately by Johnson & McKee (1971), who 
estimated T f ~ (T s /\/2) 1+ ^ as a lower limit. 4 
In this case, postshock acceleration increases the 

4 Johnson & McKee (1971) showed Tf would increase by a 



energy per particle by a factor of approximately 
r° 73 . As this is a large factor, the postshock ac- 
celeration, and the concentration of energy in the 
fastest material, are even more impressive in the 
relativistic case. 

These considerations demonstrate that the ac- 
celeration of the fluid after it is shocked is critical 
to the distribution of energy with Lorentz factor 
in the ejecta. The character of planar postshock 
acceleration is qualitatively different in the non- 
relativistic case, where the velocity increases by 
a fixed factor, from the relativistic case, where 
the final Lorentz factor becomes a power of the 
postshock value. However, in both cases there is 
a functional relationship between Tf(3f and T s f3 s . 
In reality there is a dependence on the slope of 
the density distribution as well, but this is quite 
weak and we are justified in ignoring it. Hence 
our goal is to find an approximate analytical rela- 
tionship between the final velocity of a fluid ele- 
ment and the velocity of the shock that struck it: 
one that obeys the known limits and remains suf- 
ficiently accurate in the trans-relativistic regime. 
We shall begin by examining postshock accelera- 
tion in strictly planar geometry. Then we shall 
consider how spherical expansion may reduce the 
terminal Lorentz factor of the ejecta relative to the 
planar case. 

2.2.1. Planar Geometry 

We continue with the density distribution p = 
Ph(l — r/R) n = phx n , with n = 3, which now de- 
scribes a planar slab of material of length R. We 
release shocks of varying strengths that propagate 
down this density ramp. In regions significantly 
far away from the launching region (R — r <C R), 
and except for the effects of finite resolution, the 
shock and postshock behavior ought to be the 
same among all of these runs - i.e., the function 
Tf0f(T s s ) should have the same form. This fol- 
lows from the fact that a power law density dis- 
tribution introduces no scales of its own. By ex- 
ploring the relation between shock and postshock 
velocities in simulations of different initial shock 
strength, we are able to discriminate this common 
behavior from the artifacts of the particular initial 
conditions or of the finite numerical resolution: see 



multiplicative constant if the rest mass of the fluid were 
included in the equation of state. 
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Fig. 3. — Final fluid velocity (Tf(3f) in the planar 
free coasting state compared to the velocity of the 
shock (T S (3 S ) striking a particular fluid element, 
as a function of r s /3 s . Numerical results for differ- 
ent initial energies are shown by the heavy solid 
and dotted lines. The dotted lines, at the highest 
energies, have the greatest systematic uncertain- 
ties from corrections for incomplete convergence. 
The dashed line is our expression for planar post- 
shock acceleration T f f} f /T s /3 S = C" nr + (T S (3 S )^ 
(eq. [17]). The inclined long dashed line is the 
ultrarelativistic prediction of Johnson & McKee 
(1971), T f = r^ +y5 = 0.391^+^, which we find 
underestimates Tf by a factor of ~ 2.6. 



Figure 3. 

For the purposes of numerical stability, we al- 
low the fluid to expand into a uniform medium 
of sufficiently low density that its effects on the 
accelerating ejecta (a reverse shock) can be visu- 
alized and excluded from consideration. Because 
the conversion of internal into kinetic energy is 
quite slow in planar expansion, and because of 
time-step constraints on our code, it was neces- 
sary to examine the approach of mass elements 
toward their terminal velocities in order to esti- 
mate how much acceleration would occur if they 
were allowed to expand indefinitely. A correction 
to T/3 based on the extrapolation of internal en- 
ergy to zero was made. Typical values of this cor- 
rection were < 5 — 10% for the lower energy runs 
(solid lines in Figure 3). For the most relativis- 
tic cases considered (dotted lines) the correction 
was ~ 20 — 30%. This correction procedure was 
checked and validated against the end state of a 
single simulation, allowed to expand out to a very 
large distance, at which point the internal energy 
of the zones was negligible. These corrected results 
were then further corrected for incomplete numer- 
ical convergence, for both the shock and final ve- 
locities, by considering the results as a function of 
the number of simulation zones. The largest values 
of these corrections, again for the most relativistic 
runs, were ~ 50%. 

We find that the following analytical formula 
captures the relationship between shock and ter- 
minal velocities: 

l^±^C m + (T s {3 s )^. (17) 
r sPs 

For density distributions with n ~ 3, C nr = 2.03 
(Matzner & McKee 1999). Comparing to numeri- 
cal simulations, equation (17) is accurate to ±10% 
for T S {3 S < 5, corresponding to Tf(3f < 100; and 
±25% for T s j3 s < 10, corresponding to Tfflf < 
500. The uncertainty is greatest for the most rela- 
tivistic ejecta. In the ultrarelativistic limit, equa- 
tion (17) tends towards the scaling predicted by 
Johnson & McKee (1971), but increased by a fac- 
tor of V2 1+V " ~ 2.6. 

2.2.2. Spherical geometry 

The above calculation assumes that planar 
symmetry holds throughout the postshock accel- 
eration phase of the ejecta. In fact, the planar, 
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Fig. 4. — Spherical postshock acceleration for sim- 
ulations with £i n = 0.001,0.01,0.03,0.1,0.3. Solid 
lines are regions where we have corrected for nu- 
merical convergence. Dotted lines are the highest 
resolution results without such a convergence cor- 
rection. The dashed line shows the planar result 
from Figure 3. 
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Fig. 5. — Deviation of spherical postshock acceler- 
ation from the planar limit (equation 17) as a func- 
tion of .t 1 / 3 = (1 — r/i?) 1 / 3 , for simulations with 
a power-law initial density distribution (p oc x 3 ) 
and E in = 0.001,0.01,0.03,0.1,0.3. The solid and 
dotted lines arc the same as in Figure 4. The 
long dashed line shows the nonrelativistic predic- 
tion for / sp h for an n — 3 polytrope (eq. [19]). For 
x 1 / 3 > 0.6 the chosen initial conditions no longer 
approximate a polytrope, which explains the dis- 
crepancy between theory and simulation in this 
region for nonrelativistic ejecta. Uncertainties in 
the planar solution also contribute to the accuracy 
of this prediction. 



11 



nonrclativistic flow approaches its terminal ve- 
locity extremely slowly: any fluid element must 
travel many times its initial depth before acceler- 
ation is complete. Once the fluid has traveled a 
distance of order the stellar radius, however, its 
volume is significantly larger than in the planar 
approximation, and its pressure is correspond- 
ingly less. This causes nonrelativistic acceleration 
to cease at some radius, as originally suggested 
by Litvinova & Nadezhin (1990). Comparing to 
nonrelativistic numerical simulations, Matzner & 
McKee (1999) find that the results of spherical ex- 
pansion match the results of planar acceleration, 
truncated at three stellar radii; this is also con- 
sistent with the analytical theory of Kazhdan & 
Murzina (1992). In terms of the normalized depth 
x = l- r/R, Matzner & McKee (1999) find that 
spherical truncation reduces the terminal velocity 
of a mass element originally at depth xR by the 
factor 

fsph nr (n = 1.5) = 1 - 0.51x 1/3 + 0.76x 2/3 - 1.19a; 

(18) 

/ sp h nr(n = 3) = 1 - 0.34x 1/3 + 0.24x 2/3 - 0.37a; 

(19) 

for polytropes of these indices. These results on 
postshock acceleration are valid for polytropic en- 
velopes in an outer region with rh cx = 1 — rh <C 1, 
and x 1 / 3 < 0.8. Deeper into the envelope, poly- 
tropes differ in structure if they enclose different 
core masses; this, along with changing dynamics 
deep within the explosion, causes deviations from 
these common forms (Matzner & McKee 1999). 

We have followed the postshock acceleration of 
ejecta from spherical explosions of different ener- 
gies, using for a progenitor the external power law 
distribution p = ph(l — r/R) 3 = pnx 3 , to examine 
how spherical truncation, expressed via 



sph 



spherical 
( r //3/)planar 



(20) 



changes as velocities become more and more rela- 
tivistic. Figures 4 and 5 compare the planar law 
(eq. [17]) with spherical numerical results. The 
predicted nonrelativistic truncation (eq. [19]) is 
recovered for a; 1 / 3 < 0.6, where p oc x 3 well ap- 
proximates a polytrope. At higher energies and 
velocities we observe more significant truncation in 
the outermost ejecta. Tjfif is reduced to ~ 70% of 
the predicted planar value in the most relativistic 



case. We note this estimate is somewhat uncer- 
tain, because of uncertainties in the planar solu- 
tion to which we make comparison. 

For purposes of estimating the final energy dis- 
tribution of the outer (x < 0.3, a; 1 / 3 < 0.67) ejecta 
we approximate / sp h as a constant factor, indepen- 
dent of x, varying between 0.85 and 0.7 as E- m is 
increased from 0.0 to 0.3. This behavior is cap- 
tured by 

/ sph = 0.85 - eH 2 /A, (21) 

which is approximately the value of / sp h at x = 
0.125. This expression is accurate to ±15% at any 
given x, tending to underestimate velocities at x < 
0.125, and overestimate those at x > 0.125. It has 
greater accuracy for predicting the mean value of 
/ sp h for 0.0 < x < 0.3. 

Combining expressions (14), (17), and (21), the 
final velocity of a mass element is 

T f f3 } = (0.85-^ n /2 /4Ml+P 2 ) 012 
x[2.03 + p 1 - 73 (l+p 2 ) - 21 ] 

~ (0.85 - EV 2 /4)(1.80p 5 / 8 + p 2 - 82 ) 6 / 5 ,(22) 

where the second line approximates the first 
within 1.5%. Either approximation holds within 
15% only in the outermost <~ 30% of the stellar 
radius, where equation (21) is valid. 

In an outer layer (x < 0.2), the polytropic dis- 
tribution (eq. [10]) limits to the power-law form 
(eq. [12]). The mass and density are then related 
by (Matzner & McKee 1999) 



where 



f P = 



7 P = 1 + 

n+ 1 \ w _l/(n+l) 



47T 



Ph 



(24) 
(25) 



which is Matzner & McKee's / Po (l), and we have 
defined 

« - PhR3 foa\ 
PH = -^rp (26) 

which is Matzner & McKee's parameter pi/p*. 
For n ~ 3, equation (23) is valid to within 10% 
for r > OMR and to within 50% for r > 0.75R. 
The scaling p lv oc m cx arises from the hydrostatic 
relation between mass and pressure in a thin sub- 
surface layer. 



12 



In this outer power-law region, p can be eval- 
uated by approximating r ~ R and m ~ 1 in 
equation (15): 

' - a ^(m) (27a) 

- AE^ 2 p-°- ls7 x-°- 187n (27b) 
c ^ 2 / p -°- 187 m-°- 187 ^, (27c) 

where equations (27b) and (27c) use (12) and 
(23), respectively. The error introduced by adopt- 
ing these power law approximations for p is 
roughly x, and in the outer fifth of the radius 
it is smaller than the error introduced by / sp h in 
equation (21). 

As we have noted, equation (21) for the param- 
eter / sp h, and the power law representations for 
p in equation (27), break down for x > 0.25. As 
the shock intensity p increases outward near the 
surface of the star, this implies that there exists 
a minimum value of p (hence, from eq. [22], a 
minimum final velocity) for which our equations 
accurately predict the final velocity. To estimate 
this criterion, we may set x < 0.25 and estimate 
A ~ 0.736 in equation (27b) to find 

p > 0.736 cxp(0.26n)^ i 1 n /2 / 5-° 187 . (28) 

2.3. Kinetic Energy Distribution of the 
Ejecta 

In order to estimate the observable properties of 
relativistic mass ejection from astrophysical explo- 
sions, it will be useful to know the kinetic energy 
in ejecta above some velocity E^{> Tffjf). This 
information is available through the formulae pre- 
sented in §§2.1 and 2.2. In this section, we shall 
use the above results to present more convenient 
analytical formulae for Ek(> Tfpf). We restrict 
our attention to a shell of material from the out- 
ermost layers of the star, whose mass is negligible 
(m cx = 1 — m <C 1) and through which the shock 
accelerated. We approximate the spherical trun- 
cation with eq. (21), which tends to underesti- 
mate the kinetic energy in relativistic ejecta, and 
overestimate the energy in nonrelativistic ejecta, 
although never by more than 15%. 

The fraction of the rest mass energy transmit- 
ted into kinetic energy at velocities above Ff(3f 
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Fig. 6. — Cumulative kinetic energy distribution 
of ejecta from a progenitor with p = p/j(l — r/R) 3 
for E in = 0.001,0.01,0.03,0.1,0.3. Solid lines - 
converged results from numerical simulation; dot- 
ted lines - analytic estimate (eq. 37); dashed lines 

- relativistic limit (eq. 36); and long-dashed lines 

- nonrelativistic limit (equation 35) shown for the 
two lowest energy cases. 
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IS 



Ek(>r f f3 f ) 



(T f - l)dm. 



ox: 



(29) 



where m ox refers to the ejecta traveling faster than 

r//3/. 

Lacking a relativistic theory, Matzner & McKee 
(1999) estimated the energy in relativistic ejecta as 
M re ic 2 , where M re i is the mass in which an extrap- 
olation of the nonrelativistic expressions (ignoring 
spherical truncation) predicts /3/ > 1. Examin- 
ing equations (3), (15) and (17), we find that this 
definition corresponds to p > 1/C nr : 



M rel = m ex (p > l/C nr )M e 



(30) 



Equations (14) and (17), together with allowance 
for spherical truncation (eq. [20]), show that 
the ejecta with mass M re i expand with Tff3f > 
1.18/ S ph; if / sp h — 0.85, as indicated by equation 
(21), this is Tf(3f > 1.00. M re i is thus an excellent 
estimate of the mass in relativistic ejecta. 

For a stellar progenitor with a polytropic atmo- 
sphere, Matzner & McKee (1999) find 



M i ~ in 28 - 5 - 4 - 8 /™ f-Tp ( - 
rcl - U J p \ 0.794 



10 52 erg 



2.67 7 p 



794 

1M 



5.357 P 



l-2.677p 



(31) 



where f p is defined in equation (23) and we have 
included the dependence on A; we shall present 
equivalent expressions for progenitors with radia- 
tive outer envelopes below. 

Equation (31) assumes that the shock velocity 
is a power law in depth. This assumption is valid 
in the outer ~ 25% of the stellar radius; equiva- 
lently, it holds when the value p — 1/C nr = 0.49, 
which corresponds to M re i, is above the lower limit 
in equation (28). This requires 



£ in <0.44exp(-0.52n)# 



-0.374 



(32) 



If this condition is not satisfied, the value of M re i 
in equation (31) will under predict the mass in 
material withp > 0.49, and equation (21) will over 
predict its final velocity. Even so, M ro i remains an 
accurate coefficient for higher-velocity ejecta that 
do arise in the outer quarter of the radius; c.f. 
equations (41) and (42) below. 



To evaluate M IC \c 2 as an estimate of the ki- 
netic energy in relativistic ejecta, and to deter- 
mine the distribution of kinetic energy with veloc- 
ity, we must evaluate equation (29). For nonrela- 
tivistic ejecta, equations (14), (15), (17), and (23) 
give (3f oc m 0X Qnr ^ 7p oc m ox °' 187 ^ 7p , for r ~ R and 
m ~ 1. Since E k {> T f f3 f ) = / ™ cx (3 2 f dm' e j2 for 
nonrelativistic ejecta, this implies that 



Ek(> r f (3 f ) 



2 - 0.75/ 7p 



(33) 



For extremely relativistic ejecta, the same equa- 
tions give Tf cx m cx °' 634 ^ 7p , and since Ek{> 



r f/3f) — f™ ex Tfdrh' ex in this limit, 



(34) 



Note that the energy above some final velocity 
Tf/3f is proportional to the mass m cx in ejecta 
above that velocity, quite generally. This follows 
from equations (17) and (21), which demonstrate 
that the final velocity is essentially a universal 
function of the shock velocity for a given element, 
modulo slight variations in / sp h not captured in 
equation (21). Therefore, £*(> T//3/) is propor- 
tional to mass of ejecta that was hit by a shock 
going faster than the corresponding shock veloc- 
ity r s/ 3 s . 

Eliminating m cx from equations (33) and (34) 
by means of equations (14), (17), (20) and (23), 
we find in the non- and ultrarelativistic limits 



2 - 0.75/ 7p 

/ P ^(r//3/)- (5 - 357p - 2) 

[T s d s « 1] (35) 



and 



I j-0.295 A p!/2^5.357p 

E k (>r f f ) ~ ^ '^4 

feV iHt ' l-0.63/7 P 

x /"" 7p (r//?/)" (1 ' 587p " 1) 

[T f (3 f » 1] (36) 

respectively. One formula that attains each of 
these limits is 

E k (> T ff 3 f ) = (AEll 2 f^fp*F(T f p f ), (37) 
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where 
F(T f P f ) 



2 - p f - (0.75 - 0.11/?/)^ 



(17/3/) 



-(5.35-2/ 7p )/g 



+ (r / /3 / )- (1 - 58 - 1/7p)/9 ] <?7p .(38) 

We find that setting q — 4.1 gives an excellent 
fit to the simulation results, as Figure 6 demon- 
strates. This fit is more precise than deriving 
E k {> TfPf) from equations (27c), (14), (17), and 
(21), principally because details of the spherical 
truncation are not captured in equation (21). Note 
that for progenitor density distributions that fall 
off rapidly near their edges, the energy contained 
in ejecta above a particular velocity is quite sen- 
sitive to the degree of shock and postshock accel- 
eration. For example a 10% error in the final ve- 
locity caused by equation (21) can lead to a factor 
of two error in E k {> ^fPf) in the nonrelativistic 
limit. Variations in the polytropic index n may 
cause a nr to vary by a few percent (Sakurai 1960). 
Such differences may be amplified into large vari- 
ations in E k (> T//3/), typically — 50%. Note that 
the ultrarelativistic limit is approached only very 
slowly, i.e. for > 100. This is because the 

Lorentz factor of the shock, r s , needs to be rela- 
tivistic for these limits to apply and Tff3f ^> T s f3 s 
in the relativistic regime (eq. 17). Figure 6 also 
illustrates that the cumulative energy distribution 
deviates from equation (37) in the lowest-velocity 
ejecta. This is due to the failure of several of our 
assumptions in this material: r is no longer close 
to R, / sp h decreases at lower velocities, and the dy- 
namics changes once m cx < 1 is no longer valid; 
these issues were treated in detail by Matzner & 
McKee (1999). 

Expression (37) is a remarkably simple estimate 
of the fast and relativistic ejecta from a stellar ex- 
plosion. In the following sections, we shall use it 
to predict the outcomes of explosions within ra- 
diative stellar envelopes, and of asymmetrical ex- 
plosions. We wish therefore to re-express (37) in 
a couple of forms that will facilitate these investi- 
gations. 

First, suppose one knows that the shock exceeds 
an intensity p' in a small fraction m ex (> p') 01 the 
total ejecta. What will be the kinetic energy in 
ejecta faster than Fff3f, where Tf[3f is arbitrary? 



Using equations (14), (15), and (23), we rewrite 
(37) as 

E k (> Tfpj) ~ p' 5 ' 35 ^m cx (> p')F{Y f P } ). (39) 

Here, the combination p /5 - 357p m cx (> p') is a con- 
stant independent of p' , so long as it is evaluated 
in the outermost <~ 25% of the radius (e.g., as in 
eq. 27c). This leads equations (35), (36), and (37) 
to obey the scaling (in dimensional variables) 

E k (> YfP f ) cx [iW(M cj c 2 )] 2 - 67 >M cj c 2 . (40) 

Second, we may express E k (> YfPf) in terms 
of the rest mass scale M rc \ for relativistic ejecta. 
Taking / sp h = 0.85 and p' = 1/C nr = 0.49, equa- 
tion (39) gives 



(41) 



Evaluating this at the transition from nonrelativis- 
tic to relativistic motion (17/3/ = 1), 



E k (> 1) ~ 1.23 



1 



n + 2.1 



exp 



J_\ / /sph 

3.7n/ I 0.85 



2.67 7p 



(42) 

Hence E k (> 1) ~ (0.71, 0.88, 0.94, 1.07)M rc ic 2 for 
n = (1.5,3,4,9), respectively, justifying M re \c 2 as 
an energy scale for relativistic motion. 

2.4. Envelope Structure and the Efficiency 
of Fast Ejecta Production 

An accelerating shock front produces rapid 
ejecta in much the same manner as a whip pro- 
duces its crack: by concentrating a diminishing 
fraction of its energy in an even more rapidly 
diminishing fraction of the total mass, as a dis- 
turbance travels outward. On Earth, whips are 
fashioned to yield the sharpest report for the least 
effort, so we may ask: What qualities in a stel- 
lar progenitor will maximize the energy in ejecta 
above a given velocity? 

To address this question, we rewrite equation 
(27c) with p = 0.49, as 



Mrcl 
M ei 



0.736 



2.23-Ej n 

f0.374 
Jo 



2.67~fp 



(43) 



This ratio represents the efficiency with which a 
given explosion creates relativistic ejecta. For a 
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given value of E ln , what properties of the ejected 
envelope enhance M re i/M e j? Note that ^4/0.736 
and f p are both roughly unity; moreover, equa- 
tion (27c) is validated only for E ln < 0.3. For the 
mildly relativistic explosions to which (43) applies, 
it states that M re i/M e j is a small number raised 
to the power 2.677 p . The relativistic yield of an 
explosion is therefore enhanced by making 7 P as 
small as possible, i.e., by making n as large as 
possible. For a fixed value of j p , the relativistic 
yield is also enhanced by decreasing f p , or equiv- 
alently, by decreasing the outer density coefficient 
Ph- Both of these results can be interpreted to 
say that more centrally concentrated envelopes are 
more efficient at producing relativistic ejecta. 

Second, note that the relativistic ejecta of a 
trans-relativistic explosion are enhanced - both in 
absolute terms and as a fraction of the total - if 
Ei n is increased or M e j is reduced. For a given 
progenitor mass, the relativistic yield is thus en- 
hanced by allowing as much of the star as possible 
to collapse. 

Equations (41) and (42) demonstrate that the 
kinetic energy at any final velocity Ff/3f (relativis- 
tic or otherwise) is proportional to M ro ic 2 , so long 
as this ejecta originated from a thin outer layer in 
which the shock accelerated; therefore, the above 
conclusions apply to fast ejecta even when the 
maximum ejecta velocity is nonrelativistic. 

Let us examine the criterion for relativistic 
ejecta to exist. The alternative is that the shock 
diffuses through the surface before p reaches 0.49, 
the value for which Tf/3f > 1. The optical depth 
of a nonrelativistic shock front is t s <~ 1 jp (Weaver 
1976), whereas the optical depth of the layer with 
mass M re i is r = kM ic \/(4 : ttR 2 ), where k is the 
opacity. To contain the shock front when p = 0.49, 
r > t s ~ 2, or 



R<55 



A? 



0.34 cm^- 1 10 48 erg) 



1/2 



Re- (44) 



This constraint on R, which is easily satisfied by 
the model C06 considered below {M IC \c? ~ 10 48 
erg; R <~ R & /6), becomes less restrictive as M le \ 
increases. 

The mass M re i in relativistic ejecta, and the cor- 
responding kinetic energy Ek{> 1), can be derived 
from equation (43) using the equilibrium struc- 
tures of stars with radiative outer envelopes. In 



the case of a constant (e.g., electron scattering) 
opacity, n = 3 and 

V 4 1 TIT , ■. / M \ V 2 

1 - L±/L cAd ( fj, \ j M+ ' 



f P = 0.148 



M*_ 



(At/ ^edd) 



^cdd / /Z \ 
W4 V062/ 



10 Mq 
(45) 

(Matzner & McKee 1999; note typo in their eq. 
12), where /i is the mean molecular weight in amu 
and L dd is the Eddington limiting luminosity. Us- 
ing (43) and (45), we find that the fraction of the 
ejected mass that becomes relativistic is 

3.57 / A \ 7.13 

430 M " 4 / 3 



M, 



rcl 



E\, 



M cj C 2 



A 



0.736 



JM" 1 / 3 /Vf 2 / 3 IT IT \l/3 



(46) 



(1 - i,/iedd) 4 / 3 ' 

where Mch = 1-434 Mq is the Chandrasekhar 
limiting mass for an object with /i e — 2. As 
increases relative to L c dd, the outer envelope be- 
comes significantly more dilute compared to the 
core; as we argued in the previous section, the ef- 
ficiency of relativistic ejection, M re i/M e j, increases 
correspondingly. This effect can be achieved, 
for instance, by increasing the radiative opacity 
(hence lowering L c dd): intuitively, the envelope 
must expand to carry a constant flux despite the 
higher opacity. 

Equations (42) and (46) then give 

0.88M rc ic 2 

(1.8,2.6,8.9) x 10 46 £| 2 57 L5 /3 



E k (> 1) 



M* 
6 M Q 

A 
0.736 



-4/3 



7.13 



1 - 



L c dd, 
-2.23 



-4/3 



erg, (47) 



,4 M 0/ 

using / sp h — 0.85, for stars with CO envelopes, He 
envelopes, and H envelopes of solar composition, 
respectively. In this equation, L* = 10 5 L 5 Lq and 
E in = 10 52 £ 52 erg. 

If the opacity in the outer envelope (the region 
containing M TC \) is dominated instead by bound- 
free or free-free transitions, then M le \ will ex- 
ceed the prediction of equation (47). Provided 
radiation pressure is negligible, which requires 
L* <C Add, the envelope will be a polytropc with 
n = 3.25 and 

/ „ \ 2 



f U = — 



R \ hlKoL+M*. 



GM*/j,M p 

27rfcn 
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• (48) 
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Here k — KopT 3 5 is Kramers' opacity law; we es- 
timate K ~ (37, 0.69, 1.24) x 10 24 cm 5 g~ 2 in (CO, 
He, H) envelopes (Harpaz 1994; Clayton 1983), by 
taking both relevant Gaunt factors to be ~ 0.85 
(Mrei depends only weakly on this choice.) The 
energy above T//3/ = 1 is, using / sp h = 0.85, 



E k {> 1) 



= 0.90M rc ic 2 

~ (3.2,2.4,6.3) x 10 46 £f 2 50 i5 /13 



x R 



1/13 
10 



4M B 



6 M 

2.19 



-15/13 



erg, 



A 



0.736 



6.90 



(49) 



for (CO, He, H) envelopes, respectively, where 
R = 10 10 i?io cm. In the Kramers layer, radiation 
pressure increases outward as a fraction of the to- 
tal; when it is no longer negligible, n increases 
toward seven (Chandrasekhar 1939); therefore, 
equation (49) is an underestimate in stars with ap- 
preciable L±/L c ad- In general, the larger of (47) 
and (49) gives a lower limit for the kinetic energy 
of relativistic material. 

Let us examine the criteria for equations (47) 
and (49) to correctly represent the energy in rela- 
tivistic ejecta. The first of these is given by con- 
dition (32), which states that M re i must originate 
from the outer quarter of the radius, to correctly 
give Ek(> 1). For (CO, He, H) electron-scattering 
atmospheres this becomes, using eq. (45), 



< (0.12,0.081,0.021) (° ' M) 



^odd 



V A 

Mi, 



6 M F 



4M 



0.374 

(50) 



Using eq. (48) to express the same condition 
for Kramers atmospheres with negligible radiation 
pressure, 



< 



(0.087,0.15,0.037) 



M* 

6 M CT 



0.736 
A 

15/2 



4 M G 



n 0.187 

(51) 



in the same three cases. However, recall that M Te \ 
can be used to predict the distribution of ejecta 
at sufficiently relativistic velocities (eqs. [41] and 



[42]), even when this condition indicates that it 
does not accurately determine E/-(> 1). 

In contrast, there are no relativistic ejecta at all 
if the shock breaks out of the star before p > 0.49. 
This condition, expressed above in equation (44), 
gives for Thomson atmospheres 



R 10 < (43,51,122) 



A 



1 



3.57 

0.736/ " ,J 

~ 2/3 / 6M Q \ 2/3 / 4 Mq n 1 



LcddJ 



for (CO, He, H) envelopes. In the case that 
Kramers opacity determines the structure of the 
envelope, it is still electron-scattering opacity that 
controls the width of the shock. Therefore, equa- 
tions (44) and (48) give 



R 



10 



< 



(55,48,101) 



A 



0.736 



3.64 



x L^Eir{^ 



3/5 



4M 



(53) 



2.5. Aspherical Explosions 



In equation (40), we have shown that the energy 
above some velocity scales as E^ 7lp M^ 2 ' 677p . 
Although our results so far have been limited to 
spherical or planar symmetry, this scaling sug- 
gests that asymmetries in the initial explosion may 
be exaggerated in the angular distribution of rel- 
ativistic ejecta. The dynamics of non-spherical 
explosions are more appropriately addressed us- 
ing multidimensional numerical simulations (e.g., 
Chevalier & Soker 1989; MacFadyen & Woosley 
1999; Hoflich et al. 1999; Aloy et al. 2000). Be- 
cause such simulations cannot afford the fine sur- 
face zoning necessary to capture the dynamics of 
shock acceleration, we wish to answer the follow- 
ing questions: Under what conditions may the the- 
ory developed here for spherical explosions be ap- 
plied to aspherical explosions? What is the appro- 
priate correspondence between our formulae and 
the outcome of an aspherical explosion? 

To address the first question, it suffices to note 
that, during the terminal acceleration of a nonrel- 
ativistic shock and the subsequent postshock ex- 
pansion (Sakurai 1960; Matzner & McKee 1999), 
perturbations in the postshock flow can only prop- 
agate laterally a limited distance: if the shock has 
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reached a depth x, then an acoustic wave launched 
from some point on the shock surface can only af- 
fect a region within a few times x to either side. 
We arrive at this conclusion by examining the 
self-similar solution presented by Matzner & Mc- 
Kee (1999) for the planar version of this problem. 
Equation (A2) of their appendix allows us to cal- 
culate the trajectory of a mass shell from the time 
it is hit by the shock until the final state of free ex- 
pansion. Given this solution, their equation (Al) 
gives the history of the shell's density, and hence, 
for adiabatic flow, its sound speed c s (t). The in- 
tegral of c s (t)dt gives the distance that a sound 
wave will travel along this mass shell during the 
entire period of postshock acceleration; we find 
f™^c s (t)dt = 1.17a; and 2.00a; for n = 3 and 
n = 1.5, respectively, where t n (x) is the time at 
which the shock was at the shell's initial depth x. 

This indicates that the flow will remain planar 
as long as variations in the position and strength of 
the shock front occur on length scales large com- 
pared to the current depth, so that significantly 
different patches of the surface are out of acousti- 
cal contact. This becomes even more true if the 
shock is relativistic, for then the lateral region of 
influence is reduced in extent by a factor of about 
r s by time dilation (e.g., Shapiro 1979). Like- 
wise, the onset of spherical expansion causes c s (t) 
to drop faster than the planar solution would pre- 
dict, further restricting the lateral acoustic range. 
In each direction this flatness criterion must be 
satisfied at some point as the shock approaches 
the stellar surface; after that, the motion is es- 
sentially planar. In reality, there are likely to be 
variations in the effective value of the spherical ex- 
pansion factor / sp h, if lateral expansion sets in on 
the correlation scale of the flow rather than on the 
radial scale length R; this would have the effect of 
changing the appropriate value of a; 1 / 3 . However, 
equation (21) shows that this variation can be ig- 
nored within the accuracy of our theory. 

Secondly, what is the angular distribution of 
ejecta from an aspherical explosion? Sufficiently 
close to the outside of the star that the above crite- 
rion is satisfied, each region of the surface evolves 
independently of other regions separated by more 
than a few times the initial depth. Only one pa- 
rameter controls the evolution of each region: the 
initial shock velocity. 

Suppose that a numerical simulation has pro- 



duced a shock whose intensity is p' (according to 
equation [15]) below a region with dm cx (> p')/dA. 
Equation (39) can be applied to this patch of the 
surface, provided we replace the spherical vari- 
ables Ek(> T//3/) and TO cx (r p /) with their differ- 
entials: 

— _p F(Tf/3f) , 

(54) 

this emission is in the normal direction. Note 
that here, as in (39), p' is a reference value of 
p and is not related to Tf/3f by equation (22). 
To gain insight into the implications of equation 
(54) for a given explosion, we make recourse to 
the so-called sector approximation. Here, non- 
radial motions are ignored and each angular sec- 
tor is assumed to expand as if it were part of a 
spherical point explosion with the local distribu- 
tions of mass (dM e j/dQ) and energy (dE in /dQ,). 
Equation (54) , which can be re-expressed in terms 
of dSl instead of dA, then states that the energy 
per unit solid angle above some velocity scales 
as (dEin/dn^-^idMcj/dQ) 1 - 2 - 67 ^. In the sim- 
plest case, the energy and mass are those in a sec- 
tor of a spherical explosion; in this case, the en- 
ergy per unit area (or per unit solid angle) above 
some reference velocity is unchanged. In the oppo- 
site case in which all the energy is confined within 
some small opening angle Qq, the ejecta energy 
above some reference velocity is increased by a 
factor (4tt/Q ) 2 - 67Jp . This result gives a quali- 
tative indication that aspherical explosions may 
appear far more intense (for the fortunate viewer) 
than spherical ones. Moreover, this effect is en- 
hanced in an aspherical explosion in which the 
energy injection is concentrated in directions in 
which dMcJdQ is lower than average. 

Lastly, it should be noted that Chevalier (1990) 
and Luo & Chevalier (1994) have identified an 
instability of nonrelativistic, accelerating shock 
fronts, in which the amplitude of the perturbation 
doubles every four density scale heights; it sat- 
urates in a state where density is perturbed by a 
factor 2.5 and pressure by a factor 1.3, and oblique 
weak shocks follow the leading shock. These re- 
sults were obtained for an exponential density pro- 
file. However, a shock in a power-law density 
stratification can be expected to behave similarly, 
except that the relevant perturbation wavelength 
(about two density scale heights) will shrink in 
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proportion to the distance to the surface, leaving 
behind a spectrum of artifacts in the postshock 
flow. To estimate whether the instability will satu- 
rate in a given star, note that the amplitude varies 
as p" 1 / 6 - 9 according to Luo & Chevalier (1994), 
and therefore as /3°' 78 . Because the shock's veloc- 

~ 1/2 

ity is typically p s ~ £? in when it begins to accel- 
erate (Matzner & McKee 1999), the density per- 
turbation will be enhanced by a factor ~ -E^ ' 39 
before the shock becomes relativistic (i.e., before 
p = 1). For the model of SN 1998bw described 
below, for which E ln ~ 0.003, this corresponds to 
an amplification of one order of magnitude. 

2.6. Small Ejected Masses and the Role of 
Gravity 

Thus far, we have considered only explosions 
caused by the collapse of a stellar core to a much 
denser state; this has allowed us to neglect the 
binding energy of the ejected envelope, and to use 
equations (e.g., eq. [15]) that assume a strong 
point explosion. In this section we address the 
opposite limit: mass ejection by a shock whose 
total energy is small compared to the binding en- 
ergy of the stellar envelope. This will allow us 
to analyze analytically the ejecta from collapsing 
compact objects. 

In this limit, the shock is weak as it crosses 
most of the envelope, because the thermal energy 
of a hydrostatic envelope is comparable to its bind- 
ing energy. Near the surface of the star, however, 
the temperature is proportional to the depth. A 
shock is strong if it is capable of ejecting material 
{fijC > v csc , the escape velocity), as the sound 



speed is ~ x 



1/2, 



Weaker shocks are slower 



at any reference depth; but because they become 
comparable to v csc closer to the surface, they are 
stronger (compared to the sound speed) at the 
point of ejection - so long as this occurs where 
Although we cannot use equation (15) to 
predict the shock velocity from the parameters of 
the explosion, we can nevertheless apply our for- 
mulae for shock and postshock acceleration, and 
use v esc to set the energy scale of what emerges. 

We assume that the material of interest un- 
dergoes all of its acceleration in the vicinity of 
the stellar surface; then, the ejecta undergo a 
brief phase of acceleration followed by ballistic mo- 
tion, in which they either escape (P/c > v cac ) 



or fall back. As relativistic stars are unstable, 
S g rav(env)< M env c 2 , so the explosions under con- 
sideration are necessarily nonrelativistic; likewise, 
any relativistic ejecta must escape, as v csc can only 
be a fraction of c. 

We take p to be given at some point where the 
shock has accelerated and become strong; there- 
after, the shock (eq. [14]) and postshock (eqs. [17] 
and [20]) dynamics follow the theory we have de- 
veloped. Under the assumption of ballistic motion, 
we may simply subtract the specific binding en- 
ergy, GM+I \c 2 R) = Wg SC /(2c), from the predicted 
values of T/; equivalently, the total kinetic energy 
in an external mass m cx is reduced by m ex v 2 sc /2. 
As this offset is less than the rest energy m ox c 2 , 
it clearly has no effect on ejecta with Tf 3> 1. 
However, by introducing a cutoff velocity (the es- 
cape velocity v csc ), it allows us to find the total 
kinetic energy for all the ejecta that escape. Set- 
ting m ox = Mcj and /3/c = v csc in equation (33) 
and subtracting v 2 sc M e -j2 from the result, we find 

3n GM* 

E k ~ -M ei . 55) 

5n + 8 R* 3 y ' 

Insofar as Newtonian gravity remains valid, this 
expression is accurate to within 3% even for neu- 
tron stars. 

How is the net kinetic energy in equation (55) 
distributed in velocity? For this, let us identify 
a shock intensity p csc that gives f3fC — v csc via 
equations (14), (17), (20), and (21). If we set p' = 
p csc in equation (39), then m cx (> p osc ) = 1, as p csc 
determines M e y, thus, 



E k (> r f P f ) ^p 5 cs l^F(T f P f )M ei c 2 



(56) 



which is valid so long as the binding energy may 
be neglected, i.e., Tff3f ^> v csc /c. 

For a nonrelativistic escape velocity, v csc = 
Cnr/sph?W. In equation (56), when (3f <C 1, the 
kinetic energy above some Tf/3f is independent of 

C nr and / sp h, as their influences on pl S c 5lp and 
F(Tf/3f) cancel. This is to be expected, as v osc 
sets the velocity scale and M e jV 2 sc sets the energy 
scale (eq. [55]), and shock acceleration introduces 
no scales on its own, so Ek/(M c jV 2 sc ) should be 
a pure function of {3fc/v csc . These arguments no 
longer hold for relativistic ejecta, because the ra- 
tio v csc /c becomes important, and indeed, equa- 
tion (56) depends explicitly on C nr and / sp h when 
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Tffy > 1. With C nr = 2.03 and / sph ~ 0.85, (56) 
becomes 

5.357 p 

E k (> Tffif) ~ l 43.53-9.72/„ 
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M* 6000 km\ 2 <377p 
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White dwarfs are sufficiently nonrelativistic that 
equation (57a) is a good estimate for the ener- 
getics of their ejecta. For neutron stars, the es- 
cape velocity is high enough that equation (57b) 
overestimates the kinetic energy of their ejecta by 
15 — 30%. But, the ~ 15% uncertainty in / sp h 
causes an additional uncertainty in E k {> T//3/) 
of about a factor of two, for the relativistic ejecta 
alone. 

2.6.1. Characteristic Velocities of Explosions 

Equation (55) implies that, when ejecta are 
flung ballistically from the surface of an object and 
v csc determines any fall-back, their mean velocity 
closely resembles w osc , and hence Ek ~ M e jV^ sc /2. 
Alternatively, if E m is known to exceed M e jW^ sc /2, 
for instance if E- m exceeds the binding energy of 
the entire stellar envelope, then it is not consistent 
to assume that only a small outer portion will be 
ejected. 

In general, the explosion energy is determined 
at the base of the ejecta, from which it is trans- 
ferred to the rest in a blast wave. Typically, the 
velocity scale at the base of the ejecta is the escape 
velocity of the forming remnant, or the infall veloc- 
ity of material accreting onto it, as is clearly true 
for an accretion shock that ejects material by ac- 
celerating outward. Burrows (1987) and Woosley 
& Baron (1992) identify a neutrino-driven wind, 
with a velocity <~ 0.35c (close to the escape veloc- 
ity from a newborn neutron star), as an impor- 
tant possible component of core-collapse super- 
novae and of accretion-induced collapse in white 



dwarfs; this would set an upper limit for the char- 
acteristic velocity at the base of the ejecta. In a 
thermonuclear explosion, the energy per nucleon 
prescribes a (nonrelativistic) characteristic veloc- 
ity. If the remnant is a black hole, the character- 
istic velocity may resemble that of the last stable 
orbit (mildly relativistic). However, the existence 
of relativistic radio jets and pulsar winds indicates 
that magnetically-dominated outflows violate this 
estimate (and are typically asymmetrical). But, 
the postshock pressure of a blast wave driven by 
such a wind cannot exceed the wind's ram pres- 
sure; this gives a lower characteristic velocity than 
that of the wind itself. 

One implication of this inner velocity scale is 
that equation (14) becomes invalid near the base 
of the ejecta, where it predicts higher velocities 
because it assumes a point explosion. Another im- 
plication is that the average ejecta velocity defined 
by f = 1 + E ln cannot exceed this characteristic 
inner velocity. However, the mean ejecta veloc- 
ity may be much lower, because of neutrino losses, 
photodissociation, and work done against gravity, 
or simply because a relatively massive envelope is 
swept up. 

3. Astrophysical Applications 

We now consider four astrophysical examples 
where relativistic ejecta may result from shock 
and postshock acceleration caused by an violent 
explosion in a centrally concentrated density dis- 
tribution. We use numerical simulations in tan- 
dem with our analytic approximations developed 
in §2. First we present a detailed model of an 
energetic carbon-oxygen core-collapse supernova, 
designed to model SN 1998bw, which has been as- 
sociated with GRB 980425 (Galama et al. 1998). 
The progenitor density distribution is based on a 
model kindly provided by Stan Woosley. Next we 
make a simple estimate of the ejecta properties 
of much more extreme "hypernova" explosions in 
such a progenitor, and explore their ability to ac- 
count for cosmological GRBs. Finally we present 
models of the expulsion of the outer layers of white 
dwarfs (WDs) and neutron stars (NSs), that un- 
dergo a hypothesized core collapse to a more com- 
pact state. The WD progenitor was generously 
supplied by Lee Lindblom. The density distribu- 
tions of all these models are shown in Figure 7. 
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Table 1 

Parameters of Density Models 



Model 


Ein 


A 


n 


fp 


P = PhX 3 


< 0.03 


0.68 


3 


0.63 


p = PhX 3 


0.1 


0.72t 


3 


0.63 


P = Ph^ 


0.3 


0.74t 


3 


0.63 


98bw (m ox > 2 x 1(T 5 ) 


0.003 


0.736 


9 


1.4 


98bw (m cx < 2 x 10~ 5 ) 


0.003 


0.736 


4 


0.3 


WD 


0.01-0.1 


0.705 


1.9 


0.3 


NS (m cx > 10- 4 ) 


0.01-0.1 


0.68 


1.9 


0.8 


NS (m cx < lO" 4 ) 


0.01-0.1 


0.68 


4 


0.9 



^ These values, determined by numerical simulation, arc 
~ 5 — 10% higher than predicted by eqs. (4), (5) and (6). 



3.1. SN 1998bw and GRB 980425 

SN 1998bw was a very unusual supernova. Oc- 
curring in the nearby spiral galaxy ESO 184-G82 
at z = 0.0085, equivalent to 38 Mpc for H = 
65kms _1 Mpc" 1 (Tinney et al. 1998), it is one of 
the most luminous and earliest peaking Type Ib/c 
radio supernovae ever observed (Li & Chevalier 
1999; Weiler ct al. 2000): its maximum radio lumi- 
nosity was <~ 4x 10 38 crgs s _1 , occurring 12±2 days 
after core-collapse (Kulkarni ct al. 1998). Various 
authors have derived differing shock expansion ve- 
locities from synchrotron models of this emission, 
finding (i s ~ 0.3 — 0.8 about 12 days after the start 
of the supernova (Kulkarni et al. 1998; Waxman 
& Loeb 1999; Li & Chevalier 1999). The more 
relativistic the expansion, the greater the chance 
that 7-rays could have been produced immediately 
after the star exploded. We shall consider these 
models further in §3.1.3, where we compare them 
to our predictions from the distribution of high- 
velocity ejecta. 

Because of the lack of H and He I features in its 
optical spectrum, SN 1998bw has been classified as 
a peculiar Type Ic event (Patat & Piemonte 1998). 
The optical spectra also reveal very broad emission 
line blends of Fe II, Si II, O I, and Ca II, implying 
very high expansion velocities ~ 0.2c about a week 
after the explosion. The peak optical luminosity 
(~ 1 x 10 43 ergss -1 ) is about ten times brighter 
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Fig. 7. — Comparison of initial density distribu- 
tions of our models of SN 1998bw (dotted) , white 
dwarfs (WD, dashed), and neutron stars (NS, long- 
dashed). The solid line is the idealized p = pux 3 
distribution from §2. The square marks the loca- 
tion of the outer most zone considered by Woosley 
ct al. (1999). Beyond this we extrapolate with 
polytropic models out to the circumstellar wind. 
A remnant mass of 1.78Mq has been removed from 
the center of this distribution, leaving 4.77 M Q . 
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than typical SNe Ib/Ic. Modeling SN 1998bw's 
spectra and light curve under the assumption of 
a spherically symmetric explosion of a C+O star, 
with M c j ~ (5 — 10) M Q , requires large synthesized 
Ni masses (> 0.5 M©) and explosion kinetic ener- 
gies (~ (2 - 5) x 10 52 ergs) (Iwamoto et al. 1998; 
Woosley et al. 1999; Sollerman et al. 2000; Naka- 
mura et al. 2000), again about an order of magni- 
tude greater than from typical supernovae. How- 
ever, Hoflich et al. (1999) have presented asym- 
metric explosion models that require smaller ener- 
gies. 

GRB 980425 also has some unusual properties 
distinguishing it from other bursts. Observed by 
the BATSE detector on the Compton Gamma Ray 
Observatory (Bloom et al. 1998a) and the Beppo- 
SAX satellite (Sofhtta et al. 1998), no emission was 
seen above 300 kcV, making this an example of the 
"no high energy" bursts that compose ~ 25% of 
the BATSE sample. The burst profile was smooth 
and single peaked, showing little internal variabil- 
ity. At most, only seven other bursts (0.5% of 
the total GRB population) display such profiles 
(Norris, Bonnell, & Watanabe 1999). The burst 
duration, %, was about 35 s, with the count rate 
of photons above 40 kcV peaking four or five sec- 
onds before the count rate of photons below 26 keV 
(Pian et al. 2000). Pian et al. also report that the 
supernova is coincident with an X-ray source that 
fades by a factor of two over the first six months. 

However, the most intriguing aspect of the un- 
usual events SN 1998bw and GRB 980425 is their 
spatial and temporal association (Galama et al. 
1998), occurring with < 10~ 4 chance probability 
according to a posteriori statistics. Their pecu- 
liarities only serve to strengthen the suspicion that 
they are related. Given association, the 7-ray en- 
ergy of the burst was E 7 = (8.1 ± 1.0) x 10 47 ergs, 
equivalent to ~ 10 27 g c 2 . Such a weak burst, 
thousands to millions of times fainter than the 
inferred isotropic energies of cosmological bursts, 
implies a separate GRB population. Corrobo- 
rating this suspicion, Norris, Marani, & Bonnell 
(2000) have shown that there is a power law anti- 
correlation between intrinsic luminosity and hard- 
to-soft time lags among six other bursts with 
known rcdshifts, but that GRB 980425, if asso- 
ciated with SN 1998bw, is several hundred times 
weaker than this law would predict. Moreover, 
Hogg & Fruchter (1999) have shown that if the 



rate at which GRBs occur within a galaxy is pro- 
portional to its luminosity, then the odds of detect- 
ing a burst at redshift z = 0.0085 in the current 
sample of GRBs with known redshifts are less than 
one in 10 4 . Non- Euclidean number count statistics 
limit the fraction of observed bursts resulting from 
this mechanism to < 10% (Bloom et al. 1998a). 

3.1.1. Relativistic Ejecta from SN 1998bw 

We wish to examine in detail the question of 
whether the relativistic ejecta from SN 1998bw 
were plausibly sufficient to power GRB 980425. 
Previous estimates of the mass and energy asso- 
ciated with relativistic supernova ejecta have ex- 
trapolated nonrelativistic results to the relativistic 
regime; in this manner, Matzner & McKee (1999) 
predicted sufficient energy in relativistic motion to 
account for the gamma-ray burst. Woosley et al. 
(1999) simulated shock break-out from a variety 
of progenitor models for SN 1998bw with ejected 
masses in the range 5 — 12 M and explosion en- 
ergies (4 — 28) x 10 51 ergs. Our investigation of 
SN 1998bw will be restricted to Woosley et al.'s 
model C06, which Stan Woosley has kindly pro- 
vided us. This is the bare carbon-oxygen core 
(~ 6.6 Mq) of a star that was initially more mas- 
sive (~ 25 M©); its pre-explosion luminosity is 
6.6 x 10 38 crgss -1 . The supernova forms a neu- 
tron star of 1.78 M and ejects 4.77 M with 
28 x 10 51 ergs (Table 2). It is the most energetic ex- 
plosion in the most compact, lowest mass progen- 
itor considered by Woosley et al. (1999). Among 
their models for the supernova's explosion, they 
find it gives the best fit to the observed bolomet- 
ric and spectroscopic light curves; however, signif- 
icant discrepancies, such as the fast observed rise 
in its luminosity, remain. 5 

Woosley et al. (1999) used a nonrelativistic 
code with logarithmically smooth zoning down to 
10~ 6 M . Estimating the energetics of the radia- 
tion flash associated with shock break-out, they 
concluded it was far too weak (10 42 ~ 46 erg) to 
power the GRB. Woosley et al. (1999) also es- 
timated the amount of relativistic mass ejection, 
by extrapolating their simulation to small masses 
with the shock acceleration formula of Gnatyk 

5 Iwamoto ct al. (1998) favor a more massive (M c j ~ 10 Mq) 
progenitor, with a correspondingly more energetic (up to 
~ 5 X 10 52 ergs) explosion, giving similar values for E{ n . 



22 



Table 2 

SN 1998bw Model Parameters 





6.55 M 


M rcm 


1.78 M 




4.77 M 




2.8 x 10 52 ergs 


R 


1.40 x 10 10 cm 




iog 10 r/s 

Fig. 8. — Cumulative kinetic energy distribution 
of ejecta for SN 1998bw from simulation (solid 
line) and theory: the central long dashed line is 
numerical integration of equation (29) over the 
progenitor density distribution, with T f predicted 
from equations (15), (14), (17) and (21); the up- 
per and lower long dashed lines show the effect 
of increasing and decreasing the initial energy by 
a factor of 1.5, for this method; predictions from 
equation (37) are shown by the dotted (n = 4, 
f p = 0.3, fit to outer ejecta) and dashed (n = 9, 
f p = 1.4, fit to inner ejecta) lines. The hatched 
area shows the uncertainties in the observed 7-ray 
energy of GRB 980425. 
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Fig. 9. — Density, pressure, energy (E cx ), and 
specific momentum (Tfl) of the ejecta from sim- 
ulated supernova explosion C06, after expansion 
to 9 R into a wind whose mass loss rate is M w = 
10 -5 M yr _1 . E cx is the kinetic + thermal en- 
ergy external to a particular location. The reverse 
shock due to the circumstellar interaction is vis- 
ible in the distributions of density, pressure and 
velocity. 
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(1985), given here as equation (8). On this basis, 
they estimated 10 44 — 10 45 ergs in kinetic energy 
of material with r > 3 and 10 41 - 10 42 ergs for 
r > 10. However, as trans-relativistic postshock 
acceleration (§2.2) was not included, and as we 
have improved upon Gnatyk's formula in §2, we 
wish to re-examine the topic of relativistic ejec- 
tion from model C06. 

There are four methods we may employ to ex- 
ecute this investigation. First, knowing the mass, 
radius, luminosity and composition of the progen- 
itor, we may use the analytical predictions of §2.4 
to give Ek(Tff3f > 1). For this progenitor, bound- 
free absorption dominates in an outer mass M re i. 
Equation (49), which neglects radiation pressure, 
predicts M ro i ~ 0.99 x 10 27 g with kinetic energy 
E k {> 1) ~ 8.0 x 10 47 erg. In fact as the lumi- 
nosity of C06 is about a quarter of the (electron- 
scattering) Eddington limit, radiation pressure is 
not negligible, n > 3.25 and the above estimates 
are lower limits. 

Second, if the outer density distribution of the 
progenitor were well fit by a power-law density 
ramp (eq. [23]), then with the fiducial explosion 
parameters of Table 2, we could evaluate M re i and 
Ek{> 1) directly from equations (43) and (42). 
However, for model C06 no single pair of values 
for n and f p accurately describe the whole of the 
outer region. In particular the structure changes 
around m cx ~ lx 1O~ 4 M0. The layers just interior 
to this point are well described by n ~ 9 and f p ~ 
1.4, while exterior we find n ~ 4 and f p ~ 0.3. 
To obtain the energy distribution as a function 
of velocity, 17?), we employ equation (37) 

with A = 0.736, motivated by equations (4) and 
(5), and the interior power-law fit (n = 9,f p = 
1.4). Internal variations in the density structure 
effect the accuracy of this normalization of 
Beyond the location where the structure changes 
to the n = 4 case we extrapolate using equation 
(40). This "broken power-law" estimate is shown 
by the dashed and dotted lines in Figure 8. By 
this method we find E k (> 1) ~ 2.5 x 10 48 erg and, 
from equation (43), M rcl ~ 3.0 x 10 27 g. 

A third, more accurate, method is possible: 
since we have the full density distribution of the 
progenitor model, we may directly apply equations 
(14), (17), and (21) to obtain the distribution of 
ejecta at all velocities sufficiently fast for these 
equations to hold. This technique accounts for in- 



terior density variations, which affect the nonrel- 
ativistic shock propagation and thus the normal- 
ization of the broken power-law estimate, which 
is based solely on the structure of the outer lay- 
ers. Applying this more detailed approach is com- 
plicated by the fact that the outermost zone of 
model C06 is about 2 x 10 27 g in mass; this is 
the same order of magnitude as the estimates of 
M ro i identified above. For accurate estimates, it is 
clearly necessary to continue the sub-surface den- 
sity distribution to much lower values. We there- 
fore extrapolate model C06 using a power-law at- 
mosphere with n — 4, matched to the outer edge 
of the original model. After this extrapolation, 
the smallest zone of the stellar model has a mass 
of 10 -9 M©, several thousand times smaller than 
M re i; yet, this zone is still optically thick to the 
high energy photons that drive the outward shock 
(a necessity for it to participate in shock accel- 
eration; Matzner & McKee 1999). We find that 
our results for the mildly relativistic ejecta are 
not very sensitive to our choice of n, nor to the 
exact value of the minimum zone mass. Applying 
equations (14), (17), and (21) to the augmented 
progenitor model, we find M IC \ = 3.9 x 10 27 g 
and E k (> 1) = 2.1 x 10 48 erg. These values are 
within ±30% of those obtained above by power- 
law fits. The full energy distribution is shown by 
the central long dashed line in Figure 8. The up- 
per and lower long dashed lines show the cases of 
E- m = 4.2 x 10 52 and 1.9 x 10 52 ergs, i.e. increasing 
and decreasing the fiducial explosion energy by a 
factor of 1.5. 

The final and most authoritative method to de- 
termine M re i and E k {> 1) is by direct numeri- 
cal simulation. Following Woosley et al. (1999), 
we simulate the explosion of C06 by subtract- 
ing a remnant mass of 1.78 M Q from the center 
and depositing E- m = 2.8 x 10 52 ergs in several 
pressurized, stationary zones; ignoring the effects 
of gravity, we follow the subsequent behavior un- 
til the final state of free expansion (for the rela- 
tivistic ejecta, r > 9R suffices). We embed the 
progenitor in an optically thick stellar wind de- 
scribed by the models of Springmann (1994) (see 
also Heger & Langer 1996), considering mass loss 
rates of 10~ 6 , 10~ 5 and 10~ 4 M Q yr^ 1 . The wind 
velocities are relatively low near the boundary of 
the star, increasing smoothly to achieve 90% of 
the terminal velocity, v w , within 10 — 25 R. Heger 
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& Langer (1996) set v w equal to the escape veloc- 
ity from the stellar surface, which gives a velocity 
of 3.6 x 10 8 cms -1 for our progenitor. We note 
this estimate of v w is uncertain and observed ve- 
locities from WC stars are about a factor of two 
lower (Koesterke & Hamann 1995). However, we 
find that the detailed structure of the stellar wind 
is of little consequence in our simulations of the 
explosion. After the outermost ejecta have ex- 
panded to 9 R their thermal energy has dropped 
to only a few percent of their kinetic and post- 
shock acceleration is essentially complete. At this 
point, only ~ 10 -8 M has been affected by the re- 
verse shock caused by its interaction with the wind 
(see Figure 9). This justifies the approximation 
employed by Matzner & McKee (1999), and also 
in §2 above, that the production of high-velocity 
ejecta is unaffected by the existence of circumstel- 
lar material. We shall argue in §§3.1.2 and 3.1.3 
that GRB 980425 and the radio emission from SN 
1998bw were due to the circumstellar interaction; 
however, this occurs at much larger radii, where 
the wind has reached its terminal velocity. 

Our numerical simulation places 2.3 x 10 48 ergs 
in 4.1 x 10 27 g with Tfj3f > 1, in excellent agree- 
ment with the values predicted by our shock and 
postshock theory (i.e., the third method listed 
above). The energy distribution of the simula- 
tion is shown in Figure 8 by the solid line. We 
conclude that our equations describing the shock 
speed and the postshock acceleration of trans- 
relativistic ejecta do an excellent job of predict- 
ing its energy distribution with velocity - except, 
of course, for material that has been swept into 
the reverse shock (log 10 (ry/3/) > 0.6). The upper 
and lower limits for the observed 7-ray energy of 
GRB 980425, assuming it was associated with SN 
1998bw, are shown by the hatched area in Figure 8. 
For our fiducial model, we predict ejecta traveling 
faster than T/3 > 1.35, corresponding to a velocity 
[3 > 0.8, have Ek — E 1 . Note that an explosion 
with twice our fiducial energy in a 10 M Q progen- 
itor, as considered by Iwamoto et al. (1998), has 
the same value of E in , but results in E k (> T/3) 
twice as large at any given velocity, assuming the 
same density structure. 

Compared to Woosley et al. (1999), we find two 
to three orders of magnitude more energy in ejecta 
traveling at mildly relativistic (r^ ~ 3) veloci- 
ties. This is principally because of our inclusion 



of trans-relativistic postshock acceleration, which 
boosts the final T(3 of the ejecta by a factor of 
- [2 + (T S (3 S )^]. The steeply declining density 
gradient of the progenitor, with n = 4, results in 
E k (> T f P f ) oc (r//3/)- 4 - 6 (eq. 35) in the nonrel- 
ativistic limit - a very rapid falloff with velocity. 
Increasing the final velocities of all the outer ejecta 
by even a modest factor significantly increases the 
kinetic energy above a particular ^fPf- While we 
find enough energy in mildly relativistic ejecta to 
account for the energetics of GRB 980425, the ul- 
trarelativistic ejecta arc too puny to be of con- 
sequence. We now examine if a mildly relativis- 
tic explosion can explain the observed gamma-ray 
burst and radio supernova. 

3.1.2. Early Evolution - Gamma-Ray Burst 

We assume 7-rays are produced with efficiency 
e 7 as the kinetic energy from relativistic ejecta is 
released on collision with an effectively station- 
ary circumstellar wind. For a mass m cx of ejecta 
traveling at a (mass-weighted) mean Lorentz fac- 
tor F = 1 + Ek(m cx )/ (m cx c 2 ), about half of the 
initial kinetic energy Ek{m cx ) will be liberated af- 
ter a mass ~ m ex /f of circumstellar matter has 
been swept up (e.g., Piran 1999). This estimate is 
based on energy and momentum conservation in 
the wind frame, assuming adiabatic interaction. 
Accounting for energy radiated away, leads to less 
wind material being required. For a wind of con- 
stant velocity v w and mass loss rate M w , the mass 
contained in a sphere of radius r>fl around the 
star is M w = M w r/v w . Therefore, to liberate an 
amount of kinetic energy E 1 requires an interac- 
tion with the wind material out to the radius 



r(E,) ~ 



e 7 r(r - i)m w c 2 ' 



(58) 



where e 7 cannot be much greater than a half. 

Let us suppose that the shocked shell driven 
outward by the ejecta mass m ox travels with mean 
Lorentz factor r s h c ii; this is also the Lorentz fac- 
tor of the shock that entrains the wind, if it is 
fully radiative. Ignoring the fact that this shell 
decelerates, the time for it to reach radius r is 
t' ~ r/(/3 s heiic) in the wind's frame. However, the 
observed time of the burst, t' obs 6 , is shorter by a 



l obs 



*obs/(l + 2 )> where t ohs is the actual observed time 
at Earth and z is the redshift of the burst. 
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factor 1— /Jgheii (still ignoring deceleration), so that 
r(t' ohs ) * .^ oba _ - * 2r s 2 hcll < bs . (59) 

Pshell 1 

Here, the approximation /3~^ eU - 1 ~ l/(2r 2 hell ) 
overestimates r(i obs ) for trans- relativistic motions 
(for instance, by 37% when r s h e ii = 1-7), but this 
should roughly compensate for our neglect of the 
shell's deceleration, which causes us to underesti- 
mate KCs)- 

What is the relationship between the Lorentz 
factor of the shell (L s h c ii (r j), the mean Lorentz 
factor of the ejecta that has struck the shell (L(r)), 
and the Lorentz factor of ejecta currently striking 
the shell (L/(r), say)? We must have f (r) > Tf(r) 
because slower ejecta collide with the shell at later 
times. In fact, L^(L) = Tf(Ek) is given by the 
ejecta's energy distribution, e.g., by equation (37); 
for instance, we find that to have f = 2 requires 
Tf ~ 1.7, for a variety of values of n. Also, L s h n < 
Tf, because the ejecta must be able to catch up 
with the shell. However, we expect that the latter 
quantities differ by a relatively small amount. 

The total isotropic energy of GRB 980425 in 
photons above 24 keV, assuming association with 
SN 1998bw, is E 1 = (8.1 ± 1.0) x 10 47 ergs; its 
observed duration is i bs ~ 35 s and this is ap- 
proximately equal to t' ohs , because the redshift is 
negligible. Requiring that r(t' ohs ) — r(E 7 ) inequa- 
tions (58) and (59) gives the following estimate of 
the wind's mass loss rate: 

M_ 4 3 

^1 ~ (f-l)(f/2)(r Bhc u/1.7) 2 

' (S)"(lofei)(fe)"- <60) 

where M_ 4 = M W /IQ~ 4 M Q yr" 1 and v Wt8 = 
v w /10 8 cms -1 . 

In our fiducial model for SN 1998bw, we find 
ejecta with f ~ 2 and Tf ~ 1.7 carry an en- 
ergy equal to E 7 . Ejecta with f ~ 1.6 and 
Tf ~ 1.4 contain a few xlO 48 ergs. If we take 
Tshcii ~ Tf and v w ~ 10 8 cm s^ 1 , equation (60) in- 
dicates that pre-supernova mass loss rates of order 
a fewxl0~ 4 M Q yr _1 are required to explain the 
observed 35 s duration of GRB 980425. Note that 
the radius of the circumstellar interaction is about 
400 R, far outside the region in which the wind 
undergoes its acceleration (Springmann 1994). 



Wolf-Rayet stars have mass loss rates typi- 
cally in the range 10 -6 — 10~ 4 M Q yr _1 (Hamann, 
Koesterke, & Wessolowski 1995). However, it 
is interesting that those in the the carbon-rich 
(WC) subclass are at the high end of this range. 
Koesterke & Hamann (1995) present observations 
of 25 WC stars, with a mean mass loss rate of 
~ 6 x 10~ 5 M yr _1 and about one order of 
magnitude dispersion. To produce GRB 980425 
from the mildly relativistic ejecta emerging from 
SN 1998bw, we conclude that the circumstellar 
medium be dense but within the range observed 
around Wolf-Rayet stars. 

While a dense circumstellar environment is nec- 
essary to liberate the relativistic kinetic energy of 
the ejecta, it should be checked that the 7-rays 
thus produced can escape this region to be ob- 
served at Earth. For this we calculate the pho- 
tospheric radius, R p , outside of which the 7-ray 
optical depth is r = 2/3: 

2 f°° M w h n n cs 

T = o = AnKcs / Pwdr = — — . 61 

3 J Rp 4ttv w R p 

Here, k os = 0.40/^ e cm 2 g -1 , where /i e = 
2/(1 + Xh) — 2 is the mean molecular weight 
per electron, and /k n is the relativistic (Klcin- 
Nishina) correction to the Thompson cross sec- 
tion; /kn = 0.53, 0.74, 0.84 for 7-rays of energy 
300, 100, 50keV. Thus, we have 

R p = 1.5 x 10 12 / kn ^i cm, (62) 

V w ,8 

Although this photospheric radius (~ lOOi?) is 
much bigger than the star, it is not larger than 
the interaction radius r(E 7 ) — r(t' ohs ) ~ 400i? 
identified above, implying that the GRB will in- 
deed be visible. Additionally, equation (62) pre- 
dicts smaller photospheric radii and earlier emer- 
gence for higher energy photons because of the 
energy dependence of the Klein-Nishina opacity. 
For typical conditions inside Rp these higher en- 
ergy photons are expected to be present (Weaver 
1976), and indeed the observed lightcurves of GRB 
980425 peak a few seconds earlier for higher (> 40 
keV) energies than for lower (< 26 keV) energies 
(Pian et al. 2000). 

As the emphasis of this paper is on the ener- 
getics of relativistic mass ejection, we shall not 
attempt to calculate in any detail the spectrum 
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of the observed outburst; however, note that the 
typical energies of these photons will depend on 
the mechanism of their emission. Non-thermal 
processes include inverse Compton upscattering of 
lower energy photons, or synchrotron emission by 
electrons produced in the shock front. As the lat- 
ter is the canonical mechanism for the emission 
from GRBs, we give a simple estimate of the typ- 
ical energies of these photons. Consider a shock 
of Lorentz factor r s h propagating into the wind 
(density p w = M w /{A-Kr 2 v w )). The postshock 
Lorentz factor T 2 is related to T s h by equation 
(16). The postshock energy density (Blandford & 
McKee 1976), minus the rest energy of the post- 
shock fluid, is the energy density associated with 
the postshock pressure: 

e 2 - p 2 c 2 = (4r 2 + 3)(T 2 - l) Pw c 2 . (63) 

Following Piran (1999), we assume that electrons 
and magnetic fields account for fractions e e and eb 
of this amount, respectively. The typical random 
Lorentz factor of an electron in the postshock flow 
is then 

7e = ^e.(T a - 1). (64) 

Approximating r s h c ii — T 2 , and evaluating the 
photon energy of synchrotron radiation by a typi- 
cal electron at the radius r(t' oba ) given in equation 
(59), we find 

hu syn ~ 6.3 M ?r shcll (r shcll -i) 1/2 (4r shclI + 3) 1 / 2 

With fj, e ~ 2, this gives hv syn ~ (8, 14, 29, 73) x 
(e e /0.5) 2 (e b /0.5) 1 /2(35 s/t ohs )M- 4 /v w , 8 keV, for 
r s hoii = (1-25, 1.5, 2, 3). The synchrotron model is 
thus roughly consistent with the photon energies 
observed in the burst. It suggests that the soft- 
ening of these photons over time may have been 
caused by the shell's deceleration combined with 
the declining pre-shock density. 

The earlier arrival of harder photons thus oc- 
curs naturally in both the optically thick regime, 
because R p is a decreasing function of hv, and in 
the optically thin regime, because of the decline in 
hu syn as the shell expands. 



3.1.3. Late Evolution - Radio Supernova 

Radio observations of SN 1998bw have been 
used to infer ejecta expansion velocities, assuming 
the radio emission originates in a shell just behind 
the forward shock as it propagates into the stel- 
lar wind. Kulkarni ct al. (1998) argued the lack 
of strong variability at low frequencies implies the 
shock reached a size > 10 16 cm, set by the refrac- 
tive scintillation scale of the Galactic ISM, by the 
time of the radio light curve peak at ~ 12 days. 
This implies (3 S > 0.3 averaged over this period. 
Furthermore, Kulkarni et al. (1998) pointed out 
that two constraints, one on the brightness tem- 
perature limited by inverse Compton scattering 
and the other on the total energy of a synchrotron 
emitting source, require relativistic (T/3 <~ 1.6 — 2) 
shock expansion velocities up to 60 days after the 
start of the supernova explosion. 

Waxman & Loeb (1999) and Li & Chevalier 
(1999) have presented additional models of the ra- 
dio emission. Waxman & Loeb (1999) considered 
a mono-energetic, approximately thermal, elec- 
tron energy distribution, to derive shock expan- 
sion with (3 S ~ 0.3 at t = 12 days. They argued 
the radio emission docs not imply the presence 
of a highly relativistic blast wave and that this 
weakens the link between SN 1998bw and GRB 
980425. Their model required only a weak mag- 
netic field, far from equipartition, and a high den- 
sity (~ 10 4 — 10 5 ) of emitting electrons. Li & 
Chevalier (1999) modeled a power law electron 
energy distribution and accounted for effects of 
relativistic expansion of the synchrotron source. 
They derived a range of values of the total en- 
ergy, equipartition fractions of electrons, e e , and 
magnetic fields, eb, and the ambient wind density 
that were compatible with the radio observations. 
They considered a specific, arbitrary model with 
£b = 10~ 6 and a dense (M_i/v w fi = 0.6) wind, 
which gave T S /3 S ~ 0.75 after 12 days. 

How do the above velocities compare to those 
predicted from our model? The circumstcllar 
interaction becomes nonrelativistic after a wind 
mass <~ M IC \ has been swept up; this happens at 
a radius r > 10 14 v Wi $/ M-4 cm, or after about one 
hour has elapsed since the burst. The shock is 
therefore non- or mildly relativistic (r s h e ii — 1) 
at the time of the radio observations. We assume 
the forward, radiating shock expands according to 



27 



r s oc t v , so that (3 s c — rjr/t where < rj < 1. 
The ejecta colliding with the shocked shell are in 
free expansion (/3/c = r/t), so (i s — rj(3f. The 
nonrelativistic ejecta are distributed as pf oc (3 l f 2 , 
where l p2 ~ -(8.3 + 5.3/n) (e.g., Matzner & Mc- 
Kee 1999). Because l p2 < —5, the self-similar solu- 
tion of Chevalier (1982) and Nadyozhin (1985) ap- 
plies: there, the ratio of masses of shocked ejecta 
and ambient material is a constant. Fitting the 
data in Chevalier's table, we find 



Wcx _ — Ip2 _ i o 

M w 2.7 



(66) 



within four percent. The expansion parameter rj 
is determined by the requirement that the ratio 
m cx /M w be constant: 



V 



n + 1 
L2nTT' 



(67) 



Note that the deceleration is very slow (77 ~ 1); 
this will be crucial to the persistence of rapid ex- 
pansion into the epoch of the radio observations. 

Requiring that the ambient mass within r and 
the ejecta mass traveling faster than j3f = (3 s /rj 
satisfy equation (66), and using equations (3), 
(17), (19), and (23) to specify M cj (> /?/), we find 
that for our fit (n = 4, f p = 0.3) to the outer 
envelope, 



0.43 / \ 0.13 

Vw,8 



PsheU — 0.42 [ — — 

\2.8 x 10° 2 ergs 
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■ej 



4.77 M Q 



-0.30 



1 day 



-0.13 



.(68) 



After twelve days j3 s ~ 0.30; the velocity declines 
only slowly because the amount of ejecta mass 
traveling faster than a given /3f increases steeply 
as (3f decreases. The average velocity is faster 
by only 15%, due to the slowness of the decelera- 
tion, so the mean rate of expansion is (3 S — 0.35 
at twelve days. This velocity satisfies the con- 
straints from the radio scintillation measurements, 
and is consistent with synchrotron emission mod- 
els, given their uncertainties, e.g. in the assumed 
electron distribution. Including relativistic effects 
would not much alter our result because these have 
two counteracting effects: on the one hand, they 
increase the final velocities of ejecta, thus raising 
the velocity relative to the above equation; but on 



the other, they lead to a smaller swept-up mass 
being required to slow down the ejecta, thus re- 
ducing the velocity. 

We conclude that an energetic spherical ex- 
plosion of a few xlO 52 ergs in a compact, rel- 
atively low mass carbon-oxygen core of a mas- 
sive star is capable of producing enough relativis- 
tic ejecta to account for the energetics of GRB 
980425. The results of Woosley et al. (1999) fa- 
vor this type of progenitor from fits to the light 
curve and spectra, although the explosion energy 
we require is at the high end of the range they 
considered. The timcscale of the burst requires a 
relatively dense stellar wind around the CO core, 
with implied mass loss rates <~ fewx 10~ 4 M Q yr _1 . 
Gamma-ray photons are able to escape from this 
wind to be observed at Earth. Note that opac- 
ity due to photon-photon interactions (see §3.2) 
does not constrain the minimum Lorentz factor 
of the ejecta because the observed GRB spectrum 
is very soft. These typical wind densities also al- 
low 7-rays of the correct energy to be produced by 
synchrotron emission behind the mildly relativistic 
forward shock. The shock velocity into this dense 
wind is fast enough to explain the radio properties 
of SN 1998bw days after the GRB. This simple 
spherical model is thus able to account for the ob- 
served properties of GRB 980425, if indeed it was 
associated with SN 1998bw. 

3.2. Hypernovae and Cosmological Gamma- 
Ray Bursts 

In §3.1 we found that GRB 980425 and the ra- 
dio afterglow of SN 1998bw can plausibly be at- 
tributed to the collision of mildly relativistic ejecta 
with a circumstellar wind. An interaction of this 
type has often been considered as an important 
possibility for GRB afterglows (Rees & Meszaros 
1998; Panaitescu, Meszaros, & Rees 1998; Sari 
& Meszaros 2000). External shock models have 
also been proposed (Meszaros & Rees 1993) for 
the emission from cosmological GRBs; however, 
these have been criticized (e.g., Sari & Piran 1997; 
Piran 1999) on the basis that they cannot repro- 
duce the substructure observed in many GRBs. 
This constraint does not apply to GRB 980425, 
whose profile was smooth, nor to the few other 
bursts (Norris, Bonnell, & Watanabe 1999) that 
share this property. Could cosmological bursts 
with smooth 7-ray light curves be produced by 
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Table 3 

Conservative lower limits on r in the external shock model for GRBs with afterglows 
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Yion tih iff 7 


Isotropic E-y 
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A/Tin F (on fiQ- n- — 91 


Refs. a 






(10 52 org) 


tab, =t' obs (l + z) 


[external opacity] 


[7-7 opacity] 




970228 


0.695 


0.5 


80 


9 


9 


1, 2, 3 


970508 


0.835 


0.8 


15 


17 


14 


3,4 


971214 


3.418 


30 


25 


36 


25 


5, 6 


980703 


0.966 


10 


100 


14 


12 


7, 8 


990123 


1.600 


300 


100 


28 


20 
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990510 


1.619 


30 


80 


20 


15 


3, 11 


991208 


0.706 


13 


60 


17 


14 


12, 13 


000301C 


2.03 


0.2 


10 


19 


15 


14, 15 


000418 


1.119 


4.9 


30 


15 


13 


16, 17 


000926 


2.066 


26 


25 


31 


22 
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a Table is based on data compiled by J. Greiner and S. Hempelmann at www.aip.de/~jcg/grbrsh.html. (1) Costa & Fcroci (1997); 
(2) Djorgovski, Kulkarni, & Bloom (1999); (3) Harrison et al. (1999); (4) Costa et al. (1997); (5) Hcise & in 't Zand (1997); (6) 
Mcszaros (1999); (7) from http://gammaray.msfc.nasa.gov/batsc/grb/lightcurvc/; (8) Bloom ct al. (1998b); (9) Fcroci et al. (1999); 
(10) Kulkarni et al. (1999); (11) Dadina ct al. (1999); (12) Hurley & Clinc (1999); (13) Djorgovski ct al. (1999b); (14) Smith, Hurley, 
& Cline (2000); (15) Castro et al. (2000); (16) Hurley, Cline & Mazets (2000); (17) Bloom et al. (2000a); (18) Hurley et al. (2000); (19) 
Fynbo et al. (2000). 



spherical explosions of over-energetic supernovae 
("hypcrnovae"; Woosley 1993; Paczynski 1998), in 
the same manner that GRB 980425 was produced 
by SN 1998bw? 

The cosmological GRBs with known redshifts 
require explosion energies of io 51 - 5-54 - 5 erg, i.e., 
(0.002 — 2) x M Q c 2 to be emitted as 7-rays, under 
the assumption of isotropic emission. By extract- 
ing energy efficiently from the formation of a mas- 
sive black hole at its center, a hypernova may be 
capable of such energies. (Note that the stellar- 
mass Galactic black holes are clustered around 
7 M©; Bailyn et al. 1998). However, in a spherical 
explosion of the type advocated for SN 1998bw, 
much of the energy is locked up in ejecta moving 
near the mean Lorentz factor (1 + E- m ). It is well 
known that the time scales and energies of GRBs 
require significant Lorentz factors to be optically 
thin; this would imply an even larger energy in 
low-velocity ejecta, if the necessary value of F sig- 
nificantly exceeds 1 + E- m . 

There are two potential sources of opacity that 
can shade GRB emission, and hence two con- 
straints on the Lorentz factor. The first of these 
is absorption of the observed photons in pair- 



producing interactions with other photons (Krolik 
& Pier 1991; Baring & Harding 1997); this leads to 
the "compactness problem" for GRBs. The cross 
section is highest when the center-of-momentum 
energy is somewhat above 2m e c 2 , so one must 
typically assume the number of photons varies as 
7V 7 (> v) oc v~ a , where a. > 1, to extrapolate from 
the observed photons to those that might absorb 
them (those with energies ~ T 2 (m e c 2 ) 2 / (hv) , typ- 
ically not observed). Our criterion for the fireball 
to be transparent is that r 77 < 2/3 for the ob- 
served photons. This requires that the Lorentz 
factor satisfy (Piran 1999) 

^ 2 °z"*™ s {j^){ky 2 - <«•> 

assuming the fraction of observed photon pairs 
with sufficient energy to pair produce is unity (an 
over estimate). Note, t' ohs = t bs/(l + z), where 
tobs is the observed duration at Earth and z is 
the rcdshift of the GRB. Direct observations of 
the high energy (1-20 MeV) photons from GRB 
990123 imply a ~ 2 (Briggs et al. 1999), and we 
adopt this value, as it is typical of many bursts. 
"No High Energy" bursts, such as GRB 980425, 
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Fig. 10. — Cumulative kinetic energy distributions 
of ejecta from extreme hypernova explosions, mod- 
eled in our SN 1998bw progenitor. Predictions 
are based on the analytic approximations of §2. 
The fiducial model is shown by the solid line. The 
upper and lower long- dashed lines show the effect 
of increasing and reducing E m by a factor of 1.5, 
respectively. Uncertainties in the spherical cor- 
rection factor for postshock acceleration lead to 
<~ 25% uncertainty in the final energy distribu- 
tion. The shock velocity normalization coefficient 
employed in these models is A = 0.736. How- 
ever, there is some evidence from simulation that 
this value may increase for E ln ~ 0(1). The fidu- 
cial model with A = 0.8 results in the dashed 
line. Combining these two factors, with the in- 
creased value of E- m results in a maximal model 
[dot-dashed line). The squares show the required 
isotropic energies E 1 /e 1 , with e 7 = 0.5, of the 
cosmological GRBs with detected optical after- 
glows, plotted at the minimum values oiTfPf im- 
plied by external (solid) and 7 — 7 (open) opac- 
ity constraints (Table 3). Note that the minimum 
Tf ~ r/2. The dotted arrow shows the effect on 
the 7 — 7 opacity constraint for GRB 990123 of 
adopting a smaller E 1 = 1.2 x 10 54 ergs (e.g. Briggs 
et al. 1999), consistent within cosmological uncer- 
tainties, while the dashed arrow shows the effect 
of applying the constraint to the first major pulse 
of the lightcurve with FWHM ~ 10 s and energy 
- 1/2 of the total. 



with a > 4.5 between 100 and 300 keV (Pendleton 
et al. 1997), are not severely constrained by this 
opacity. There is some evidence for harder spectra 
(a <~ 1 — 1.5), up to 100 MeV - GeV energies in a 
few bursts where observations are available (over 
the same time periods as the keV burst) (Hurley et 
al. 1994; Schneid et al. 1995; Dingus et al. 1997), 
and such hard slopes would increase the minimum 
Lorentz factor. The value of T in equation (69) 
ought to represent a frame in which the photon 
distribution is isotropic, i.e., the frame of the cool- 
ing postshock fluid; for simplicity, we shall identify 
this with the mean (mass- weighted) Lorentz factor 
of the relevant portion of the ejecta, f . 

In an external-shock model, the surrounding 
material whose presence is necessary to liberate 
the burst energy inevitably introduces a second 
source of opacity. Adopting a wind model (p oc 
r~ 2 ) for this material, and requiring its optical 
depth be less than 2/3 (i.e., requiring r(E 1 ) = 
r(t' obs ) > R p in equations [58], [59], and [62]), we 
find 



f (f - l)lt hell > 



3-Ky/knK CS 



Note, Tf ~ r shcll ~ f(> r,)/2, for n 
Tf » 1, (equation 34), and so 



r > 11 
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10 52 



erg e 



1/6 
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-1/3 



(70) 



4 and 



(71) 



as long as the implied r> 1. 

Often, models of GRBs or their afterglows 
assume a uniform rather than a wind ambient 
medium. In this case, the opacity of the swept-up 
gas is not at issue, if the burst is not in an ob- 
scured region. In an external shock model, one 
may constrain f by assuming a hydrogen number 
density nu and requiring that sufficient mass be 
swept up to thermalize the kinetic energy of the 
ejecta: 



E*y 



0.5 



10 52 



erg e n 



1/8 



'•obs 

35s 



-3/8 



-1/8 



(72) 

This constraint on T is less severe in more dense 
regions, which is why the opacity constraint is the 
more relevant one for wind interaction models. 

In an external-shock model, t' ohs is the total du- 
ration of the burst, and significant substructure 
is difficult to explain (Sari & Piran 1997); in an 
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internal-shock model, t' ohs is the duration of an in- 
dividual pulse. In either case, the total duration 
sets a lower limit for the Lorentz factor: we give 
these lower limits in Table 3, for ten bursts with 
known rcdshifts. Like GRB 980425, most of these 
bursts show little fine (sub-second) structure in 
their 7-ray light curves. The profiles are quite sim- 
ple with most of the energy contained in one or two 
pulses lasting > 10 s. Note that some substructure 
may be produced by the interaction of an external 
shock with an inhomogeneous medium (Dermcr 
& Mitman 1999; Fenimore et al. 1999). Instabili- 
ties in accelerating shocks (Luo & Chevalier 1994) 
may also create variability. Strong density jumps 
in the progenitor may lead to variable energy re- 
lease from the postshock flow as different layers 
crash into the decelerating region. Structure was 
observed in the radio emission from SN 1998bw 
(Kulkarni et al. 1998; Wieringa, Kulkarni, & Frail 

1999) , which is indubitably an external shock; the 
same can be said for fluctuations in other radio 
supernova light curves (e.g., 1979c; Weiler et al. 

2000) , which are attributed to variations in the 
ambient density. 

To gauge how difficult it is for a spherical hy- 
pernova to satisfy the above constraints on energy 
and velocity of relativistic ejecta, we adopt an ex- 
treme explosion of E- ln = 5 x 10 54 erg, perhaps 
liberated during the formation of a 10 M black 
hole. We note that models for magnetic extrac- 
tion of this energy on timescales short compared 
to the break-out time may require field strengths 
in excess of 10 15 G (MacDonald et al. 1986). For 
lack of a better model, we surround the nascent 
black hole with the ~ 5M Q envelope of model C06 
(§3.1). This gives E ln ~ 0.59, placing the explo- 
sion in a regime untested by our simulations. Re- 
gardless, we employ our analytic approximations 
to investigate the expected energy distribution of 
the hypernova's ejecta (Figure 10). 

We compare these predictions with the required 
ejecta energies and minimum Lorentz factors from 
opacity constraints of GRBs with observed op- 
tical afterglows, listed in Table 3. Note that 
Tf ~ f (> Tf)/2, for n — 4 and Tf » 1, (equation 
34). We require E k (> f/2) = E 1 /e 1 . There is 
a wide range of observed energies, yet quite simi- 
lar minimum Lorentz factors are implied from the 
photon-photon and external opacity constraints, 
spanning a range from about 10 to 40. These 



values are lower than those commonly quoted in 
the literature (e.g. Lithwick & Sari 2001), because 
the bursts in question have relatively smooth time 
profiles, and because our adoption of the external- 
shock hypothesis dictates that the age of a burst, 
rather than the timescale of its substructure, be 
used to constrain the Lorentz factor. Our extreme 
hypernova models can satisfy these requirements 
in three, and potentially six (given the uncertain- 
ties), of ten cases. However, the most energetic 
bursts require preferential viewing of asymmet- 
ric explosions (§2.5). Since the most luminous 
GRBs do not differ from typical GRBs in any 
significant manner except luminosity, it is there- 
fore likely that asymmetric explosions are required 
for all cosmological GRBs. Nonetheless, these re- 
sults demonstrate the potential importance of the 
trans-relativistic shock acceleration mechanism for 
GRB engines. 

Afterglows potentially pose a difficult con- 
straint on spherical models invoking shock ac- 
celeration and external emission, as there must 
be much more energy in slower ejecta than im- 
plied by the burst. We have found Ek(> Tff3f) cx 
l/(Tf/3f) (equation 36), approximately, for values 
of TfPf significantly above 1 + E ln ; the energy 
rises even more sharply for low values of Tf(5f 
(Figure 6). As the emitting region is continuously 
"refreshed" as it decelerates by this increasing 
source of energy (Rees & Meszaros 1998; Sari & 
Meszaros 2000), afterglows should evince these 
even higher energies. This is especially true if the 
emitting shell makes a transition from radiative to 
adiabatic dynamics between the burst and its af- 
terglow. Frecdman & Waxman (1999) argue that 
afterglow observations indicate an energy com- 
parable (within a factor of three) to the energy 
emitted in 7-rays, despite an assumed decelera- 
tion from r > 100 to T ~ 10. Moreover, the 
emitted 7-ray energy places a lower limit on the 
energy content (per solid angle) of the emitting 
gas in its early phase. 

We use the results of Sari & Meszaros (2000) to 
evaluate the severity of such afterglow constraints. 
Assuming a wind ambient medium p oc r~ 2 and ul- 
trarelativistic ejecta from an envelope with poly- 
tropic index n, equation (34) and equation (3) of 
Sari & Meszaros gives a swept-up kinetic energy 
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that varies as 

^«4 2 s 7+n)/(2 ' 7+7 - 7 "\ (73) 

i.e., E k oc ^ b ° s 2 ' - 17 ^- 13 ) for „ = (4,7,oo). The 
increase in kinetic energy from 35 seconds to ten 
days is thus a factor of ~ (7.6, 5.6, 3.7) in a model 
involving a hypernova and its wind, for these three 
values of n. Considering the uncertainties in the 
Freedman & Waxman results, we conclude that it 
cannot yet be argued that afterglow calorimetry 
rules out approximately spherical hypernovae as 
the origins of some cosmological GRBs, especially 
those with relatively low E 1 for which Freedman 
& Waxman find the afterglow to be about three 
times more energetic than the GRB. 

All of the above constraints are somewhat as- 
suaged if the progenitor or its explosion are asym- 
metric. As discussed in §2.5, an assessment of 
the implications of asymmetrical explosions re- 
quires that multidimensional simulations be inter- 
preted using our equation (54). Under the admit- 
tedly simplistic sector approximation, we found 
that relatively small variations in the energy per 
unit mass of the ejecta can have substantial ef- 
fects on the energy of the resulting relativistic 
ejecta. More accurately, given a shock intensity 
p below a given column density on some patch of 
the stellar surface, the kinetic energy in that direc- 
tion (above some chosen velocity) varies asp 5 ' 357p . 
For n = (4, 7), the shock intensity p need only be 
enhanced a factor (1.99,2.12) at some reference 
column to increase the kinetic energy in that di- 
rection by a factor of one hundred, at all veloci- 
ties for which the ejecta experienced an accelerat- 
ing shock. This enhancement would be enough to 
bring even GRB 990123 within the reach of hyper- 
nova models. We eagerly await simulations that 
will evaluate whether this degree of asymmetry is 
realistic. 

In summary, while our ad hoc and extreme 
spherical hypernova model is consistent with many 
of the GRBs with observed afterglows, it appears 
likely that the most energetic bursts require asym- 
metric progenitors or energy injection. Given this 
possibility, trans-relativistic blast waves from hy- 
pernovae may in principle account for the ener- 
getics and ejecta velocities of cosmological bursts, 
though it is too early to decide if this mechanism 
is primarily responsible for their origin. The as- 
sociation of these cosmic GRBs with the stellar 



population of their host galaxies (Bloom et al. 
2000b), and more tentatively with star forming 
regions (Kulkarni et al. 2000; Djorgovski, Bloom, 
& Kulkarni 2000; Galama & Wijers 2000), argues 
in favor of a connection with the deaths of mas- 
sive stars and against scenarios involving neutron 
star mergers. The difficulty of avoiding baryon 
loading while trying to escape the stellar progen- 
itor, encountered even in highly asymmetric jet 
models (Aloy et al. 2000), is overcome in the hy- 
drodynamic shock acceleration mechanism. The 
efficiency of 7-ray production is higher for a blast 
wave producing these photons from both external 
and internal shocks, compared to one that radiates 
solely via internal shocks (Kumar 1999). While 
external shocks are inefficient at producing sub- 
structure (Sari & Piran 1997), the lightcurves of 
the bursts with observed afterglows are relatively 
smooth, with little power in sub-second features. 
It remains to be demonstrated if our model can ac- 
count for moderate variability, perhaps via the in- 
teraction with circumstellar inhomogcneitics (Dor- 
mer & Mitman 1999; Fenimore et al. 1999) or via 
the production of postshock velocity perturbations 
by the growth of instabilities in the accelerating 
shock (Luo & Chevalier 1994) or by the pres- 
ence of steep density shelves in the progenitor's 
atmosphere. Further investigation is also required 
to determine whether a model with monotonically 
decelerating forward and reverse shocks produces 
spectra (see Sari & Meszaros 2000) that resem- 
ble those of cosmological GRBs as well as their 
afterglows. 

3.3. Ejecta from Compact Objects 

3.3.1. White Dwarfs 

Accreting white dwarfs near the Chandrasekhar 
mass face two quite distinct fates. Type la super- 
novae involve the release of ~ 10 51 ergs in 1.4 M Q 
(.Bin ~ 4 x 10~ 4 ), in either a deflagration or a 
detonation (Woosley & Weaver 1986). Approx- 
imating a detonation as an adiabatic point ex- 
plosion, equation (37) with f p = 0.3, n = 1.9, 
and A = 0.705 (values derived from a model pro- 
vided by Lee Lindblom) gives Ek(Tf > 3,10) = 
(35, 1.1) x 10 39 ergs, respectively (see also Berezin- 
sky et al. 1996, who presented an approximate 
treatment of shock and postshock acceleration in 
WDs). These energies are too small to produce 
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Fig. 11. — Kinetic energy distribution with ve- 
locity for the outer 0.1 M of ejecta from white 
dwarf explosions. Numerical simulation results are 
shown by the solid lines, and analytic predictions 
from §2, integrated over the progenitor density dis- 
tribution, by the long dashed lines. The agreement 
is excellent for the outermost ejecta, but, as ex- 
pected, deviates towards the interior, where ap- 
proximations of the analytic theory break down. 
Models of AIC (Fryer et al. 1999) tend to place 
- 10 50 ergs in 0.01 - 0.1 M . 



observable GRBs. 

Another possible outcome, if the central den- 
sity is high enough to allow electron capture prior 
to carbon or neon ignition, is accretion-induced 
collapse (AIC) from white dwarf to neutron star 
(e.g., Woosley & Baron 1992; Fryer et al. 1999, 
and references therein). Such events resemble the 
core-collapse supernovae of massive stars, but lack 
the ram pressure and long timescale introduced 
by an extended infalling envelope. Goodman et 
al. (1987) suggested that neutrino annihilations 
above the neutrinosphere of a newly-formed neu- 
tron star could vent <~ 0.3% of the neutrino lumi- 
nosity into the nearby volume; if this region were 
baryon-free, a relativistic wind of energy ~ 10 50 
erg would result, which could drive a GRB (Dar 
et al. 1992). However, the neutrinosphere under- 
lies a massive layer of baryons (Woosley & Baron 
1992), so there results instead a neutrino-driven 
wind of velocity <~ c/3. This may cither expel 
an exterior envelope, or expand freely if there is 
no envelope left (Fryer et al. 1999). There is not 
sufficient energy in such winds to power the cos- 
mological GRBs listed in Table 3, and moreover 
they cannot satisfy the Lorentz factor constraint, 
(69) or (71). They could only produce fast ejecta 
through shock acceleration; but, since the mean 
velocity of the ejecta cannot exceed c/3 (§2.6.1), 
only a small fraction of the energy is available at 
appreciable Tf. 

In addition to the energy injection from this 
wind, there is also a prompt shock from the core 
collapse (which stalls; Fryer et al. 1999), and "de- 
layed" heating by neutrino absorption behind the 
stalled shock; the latter drives ejection in the 
model of Fryer et al.. This simulation has ~ 10 50 
ergs deposited in ~ 0.1 M Q , for a mean expan- 
sion velocity of ~ c/30. Equation (55) would pre- 
dict this total energy for R <~ 1500 km. In the 
calculation, (20%, 60%, 90%) of the ejected mass 
fell below (200, 300, 400) km before being ejected 
(C. Fryer, 2000, personal communication). This 
indicates that the theory of §2.6, which neglects 
infall motions and assumes ejection from a thin 
outer layer, is not directly comparable to these 
AIC models. 

Because the total energies are too low and the 
typical expansion velocities too slow, Woosley & 
Baron (1992) and Fryer et al. (1999) conclude that 
AICs do not cause cosmological GRBs like those 
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listed in Table 3. Although this is a valid con- 
clusion, for completeness we investigate the rela- 
tivists ejecta from these events, including cases 
of more extreme energy release than predicted by 
Fryer et al. (1999). 

Our 1.38 M Q progenitor model, from Lindblom 
(1999), uses for its equation of state a fully degen- 
erate, non-interacting Fermi gas. The initial den- 
sity distribution is shown in Figure 7. We note 
this model has a radius twice that of the WD con- 
sidered by Fryer et al. (1999). 

Figure 11 shows E k (> T/3) vs F/3 for explo- 
sions of different energies in our WD progenitor. 
As these simulations neglect gravity and the ini- 
tial collapse motion of the progenitor, they should 
be considered rough approximations to the re- 
sult. The lowest curve in this figure matches 
the mass and energy of the Fryer et al. simu- 
lations; higher curves address hypothetical cases 
with much higher energies. Although an asym- 
metrical explosion, which may be expected in the 
collapse of a rapidly rotating object, might en- 
hance the energy along the line of sight, we con- 
cur with Woosley & Baron (1992) and Fryer et al. 
(1999) that AICs do not make cosmological GRBs. 
We are not in a position to judge whether magne- 
tized winds or jets from such events (Shaviv & Dar 
1995) produce GRBs. 

3.3.2. Neutron Stars 

Explosions in neutron stars (NSs) powered by 
a phase transition to a "strange" star have been 
considered as a potential GRB mechanism (Fryer 
& Woosley 1998, and references therein). Fryer & 
Woosley simulated the collapse and bounce result- 
ing from a neutron star phase transition to derive 
properties of the ejecta. These were then used to 
implement higher resolution simulations of shock 
propagation in the outer 0.01-0.001 M Q , with reso- 
lution down to 10 -11 M©. This latter stage made 
use of a nonrelativistic code, with a simple ap- 
proximation based on energy conservation to de- 
rive upper limits on the kinetic energy distribution 
at high (> 40) Lorentz factors. However, this ap- 
proximation neglects the concentration of energy 
in the outermost ejecta at the expense of the inner, 
due to pressure gradients in the flow (see §2.2). 

Fryer & Woosley's most energetic explosion in- 
volved 0.017 M Q ejected with ~5x 10 51 ergs from 
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Fig. 12. — Kinetic energy distribution with veloc- 
ity for ejecta from neutron star explosions, after 
subtracting off the initial gravitational binding en- 
ergy. Results for two simulations are shown (solid 
lines), and compared to analytic predictions from 
§2, integrated over the progenitor density distri- 
bution (long dashed lines). The analytic predic- 
tions amplify the dependence of velocities result- 
ing from abrupt changes in the initial density dis- 
tribution compared to simulation, which can more 
accurately handle these effects. The analytic pre- 
dictions are used to further explore the parameter 
space. Model C, in which ejecta of mass 0.017 M Q 
carry away 4 x 10 51 ergs of kinetic energy (solid 
square) is closest to the energetic model of Fryer 
& Woosley (1998). The solid triangles are Fryer 
& Woosley's predicted upper limits from nonrel- 
ativistic shock propagation models, including a 
relativistic correction. Our estimates, including 
the direct energy approximation of equation (57b) 
(open triangles) are about two orders of magni- 
tude higher. The dotted lines show mass contours. 
These reveal the greater ejecta masses that result 
from the explosions of higher energy. The open 
squares show the 7 — 7 opacity constraint of the 
three weakest cosmological GRBs from Table 3. 
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a neutron star of initial mass 1.4 M . Using the 
Gnatyk shock formula (eq. 8) and their simple es- 
timate of postshock acceleration, they inferred an 
upper limit of 10 46 ergs for material traveling with 
Tf > 40, and concluded that NS phase transi- 
tions cannot power cosmological GRBs. However, 
as these events are sufficiently energetic to explain 
some of these bursts, and as the mean velocity im- 
plied by their ejecta is ~ 0.51c (i.e., f ~ 1.16), we 
wish to revisit this issue with our improved treat- 
ment of the dynamics. 

As the injected energy is small compared to the 
binding energy of the NS envelope, and as only a 
small fraction of this envelope is blown away, the 
results of §2.6 apply. In the simulation of Fryer 
& Woosley (1998), the stellar surface is shocked 
after it has fallen from 11 km to 10 km. We use 
the low-density Harrison- Wheeler (HW) equation 
of state (Glendenning 1997) to obtain the appro- 
priate value of n underneath 0.017 M Q of mate- 
rial: the pressure at this depth is 5 x 10 32 dyne 
cm~ 2 , so we adopt the value n ~ 12 appropriate 
for the highest pressure listed (10 31 dyne cm -2 ). 
Using these parameters, equation (55) predicts a 
hnal kinetic energy of 3.3 x 10 51 erg, sufficiently 
close to the value (<~ 5 x 10 51 erg) quoted by Fryer 
& Woosley (1998) to give us faith in our approx- 
imations. As an added check, note that Fryer 
& Woosley consider a weak explosion that ejects 
1O~ 3 M with - 3 x 10 50 erg; equation (55) pre- 
dicts 2.0 x 10 50 erg for this ejecta mass. Thus, we 
see that the theory of §2.6 underestimates the ki- 
netic energy in the simulations by ~ 50%. This is 
not due to our approximation of Newtonian grav- 
ity as Fryer & Woosley's calculations were also 
Newtonian, but may be due to our neglect of in- 
fall of the outer layers, before they are shocked, or 
their neglect of special relativity. 

We apply equation (57b) to predict the kinetic 
energy in ejecta at some higher Lorentz factor. 
Note a value of n = 4, appropriate for the outer 
distribution, is used in evaluating F(Tf/3f), while 
n = 12, appropriate for the base of the ejecta, is 
used to evaluate the preceding terms that normal- 
ize the distribution. Taking / sp h = 0.85, we obtain 
kinetic energies of (8.2,2.4,0.85,0.45,0.35) x 10 49 
ergs for Tf > (5, 10,20,30,40); presumably, these 
energies should be increased by ~ 50% to account 
for the discrepancy noted above. This gives an 
estimate of <~ 0.52 x 10 49 erg in material with 



Tf > 40, more than two orders of magnitude 
greater than Fryer & Woosley's inferred upper 
limit. 

We validate this result with direct simulation of 
relativistic explosions in a neutron star progenitor, 
constructed using the Harrison- Wheeler equation 
of state at low densities (Glendenning 1997) and 
with the pressure extrapolated to higher densities 
with a polytrope, using Newtonian gravity. The 
central pressure and polytropic index were var- 
ied to reproduce a realistic mass and radius. The 
initial density distribution is shown in Figure 7. 
Our method follows that employed for explosions 
in white dwarfs, except now it is necessary to ac- 
count for the progenitor's gravity (§2.6). We make 
a first order correction to our non-gravitating sim- 
ulations by subtracting off the initial gravitational 
binding energy of ejecta from their final energy 
distribution. In Figure 12 we compare simulation 
results with our analytic approximations to the 
shock dynamics, applied to the progenitor den- 
sity distribution, again accounting for the gravi- 
tational binding energy. We also show (open tri- 
angles) the direct energy approximation (equation 
57b), without the calibration correction to Fryer 
& Woosley results, which acts to increase the en- 
ergies. Our different methods are in good agree- 
ment, and are larger than Fryer & Woosley's es- 
timate (solid triangles) by between two and three 
orders of magnitude, a fact which highlights the 
importance of a correct treatment of postshock ac- 
celeration. 

For a uniform ambient density of density nu ~ 
1 cm~ 3 , equation (72) indicates that Tf ~ T/2 > 
44 if a GRB with E 1 = 10 52 ergs is to last ~ 
35 s in its rest frame, although lower values of 
Tf may be acceptable if the burst occurs in the 
dense circumstellar environment of a companion. 
The photon-photon opacity constraint, eq. (69), 
requires Tf ~ t/2 > 5 for such a burst, if the 
photon index is a = 2. 

Considering the cosmological GRBs listed in 
Table 3, we see that neutron-star phase transi- 
tions can match the energies of several such bursts 
(970778, 970508, and 000301C) if Lorentz factor 
constraints are ignored. However, the required 
Lorentz factors seem to be about > 5 — 10 for 
the slowest part of the ejecta, if in a sufficiently 
dense environment. The kinetic energy of such 
ejecta is only a few percent of the total, about 
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10 50 ergs. These considerations argue against NS 
phase transitions as the origins of cosmological 
GRBs, unless the ejected mass has been signifi- 
cantly underestimated or the energy can be colli- 
mated. 

4. Conclusions 

Any viable model for the central engines of 
gamma-ray bursts, especially for those bursts 
known to be at cosmological distances, must ex- 
plain the existence of large energies (10 51-54 erg, 
if isotropic) in ejecta moving sufficiently fast that, 
when converted into radiation, it escapes as non- 
thermal radiation. Such velocities require very lit- 
tle "baryon loading," a constraint that many the- 
ories for gamma-ray burst central engines address. 
In this paper, we have investigated the possibil- 
ity that shock hydrodynamical acceleration, first 
proposed as a source of GRBs by Colgate (1974), 
may concentrate sufficient energy in the fastest 
ejecta to satisfy this constraint. We consider or- 
dinary and over-energetic supernova explosions, 
making specific models for the hypothesized class 
of supernova gamma-ray bursts (S-GRBs; Bloom 
et al. 1998a) whose existence was suggested by 
SN 1998bw and GRB 980425. We also calculate 
the yield of relativistic ejecta from hypothetical 
models for extraordinary "hypernova" explosions, 
as well as possible accretion-induced collapses of 
white dwarfs and sudden condensations of neutron 
stars. 

As an exemplary case, we have considered in 
detail (§3.1) specific models for supernova 1998bw 
and the gamma-ray burst (GRB 980425) that it 
may have produced, considering the constraints of 
the GRB and the supernova's early radio emission 
in the context of progenitor models constructed 
(by Woosley et al. 1999) to fit its light curve. One 
goal of this investigation has been to evaluate as- 
sertions (Wang & Wheeler 1998; Woosley et al. 
1999; Hoflich et al. 1999) that spherical models 
for this explosion cannot explain the GRB. After 
properly accounting for the the trans-relativistic 
shock and postshock dynamics as well as the inter- 
action with a surrounding stellar wind, we are able 
to account for all of the observations of this event 
with a spherical energetic Type Ic supernova from 
a progenitor with high (~ few x 10~ 4 M Q yr _1 ) 
mass loss rate, consistent with observed carbon- 



rich (WC) Wolf-Rayet stars. Our conclusion that 
the fastest ejecta need not be asymmetric is remi- 
niscent of the fact that the radio shell around SN 
1993J was spherical (Marcaide et al. 1997) despite 
the assertion (Hoflich et al. 1996) that its polar- 
ization implied an asymmetrical explosion. Our 
model for SN 1998bw and GRB 980425 invokes a 
circumstellar environment similar to that inferred 
around SN 1997cy, which may be associated with 
GRB 970514 (Turatto et al. 2000); however the 
two-second duration of this burst may be difficult 
to reconcile with the external shock hypothesis. 

In order to quantify the process of relativis- 
tic mass ejection, we have extended the nonrel- 
ativistic theory of supernova explosions (Matzncr 
& McKee 1999) into the relativistic regime. Find- 
ing that the energetics of such ejecta depend sen- 
sitively on the velocity of the shock front as it 
emerges from the stellar surface, we present in §2.1 
a formula for this velocity (eqs. [14] and [15]) that 
handles the transition from nonrelativistic to rel- 
ativistic flow more accurately than does the pre- 
diction of Gnatyk (1985) (equation 8). Whereas 
for Matzncr & McKee (1999) it sufficed to take a 
universal value for the normalization of the shock 
velocity (and also applied an overall energy conser- 
vation constraint to their ejecta-density models), 
we have found it necessary to present (eqs. [4], [5], 
and [6] ) a refined estimate for this coefficient - one 
that responds to the initial density distribution of 
the progenitor. 

The motivation for such precision is the sensi- 
tive dependence of the final velocity of a mass ele- 
ment on the velocity of the shock that struck it: in 
the relativistic, planar limit this is T f cx T^' 7 . This 
postshock acceleration, neglected in many works 
that invoke shock acceleration, is crucial for the 
production of appreciable energy in fast ejecta. In 
§2.2 we present formulae for this postshock accel- 
eration in the planar limit (eq. [17]) and for the 
effect of spherical geometry to reduce the acceler- 
ation (eqs. [20] and [21]). The latter is somewhat 
uncertain, due to the limitations of our numeri- 
cal calculations; however, since this correction is 
a factor of order unity, the energetics of relativis- 
tic ejecta are not especially sensitive to it. Be- 
cause we account properly for postshock acceler- 
ation, our predictions for the energetics of rela- 
tivistic ejecta are much larger than presented in 
previous works that considered hydrodynamical 
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shock acceleration in the context of supernovae, 
accretion-induced collapses of white dwarfs, and 
phase transitions in neutron stars. 

With the dynamics of trans-relativistic mass 
ejection well in hand, we turned in §2.3 to an eval- 
uation of the amount of kinetic energy expected 
in ejecta traveling above some given lower limit. 
Equation (37) represents an accurate approxima- 
tion for this kinetic energy distribution under the 
assumption that the outer stellar density distribu- 
tion is that of a polytropic atmosphere of index 
n. This formula shows that the kinetic energy 
distribution rolls over quite slowly (fig. 6) from 

E k OC ^(3-36+5.36/^ tQ ^ K r -(0.68 + 1.68/n) j &g 

the final velocity in question goes from nonrela- 
tivistic {T f (3f < 1) to relativistic (T f /3 f > 1); this 
slow limiting behavior comes from fact that rela- 
tivistic final velocities originate from only mildly 
relativistic shocks. For typical values of n, the 
kinetic energy in relativistic ejecta decreases as 
~ 1/r 1 /; this gentle decline with velocity allows 
for significant energy at relatively high T / . 

We next developed (§2.4) a theory to assess 
what characteristics of a stellar envelope make it 
especially efficient at converting a given dose of 
injected energy into relativistic motion - the as- 
trophysical analog of how to design a good bull- 
whip. If the envelope mass is fixed, we found that 
more centrally-condensed envelopes arc more effi- 
cient producers of high-velocity ejecta. If instead 
the progenitor star is given, then the energy in 
fast ejecta is enhanced by allowing as much of it 
as possible to collapse (minimizing M c j), as this 
increases the energy per mass of what escapes. 

Evaluating the energetics of relativistic ejecta 
for radiative stellar progenitors in §2.4, we derived 
expressions appropriate for regions dominated by 
electron scattering (eq. [47]) or bound- free or free- 
free (eq. [49]) absorption, assuming different pos- 
sible compositions for the outer envelope. For the 
Thomson opacity case, the central concentration 
of the atmosphere - and therefore its efficacy at 
producing relativistic ejecta - is enhanced if the 
stellar luminosity approaches the Eddington limit. 

It should be kept in mind that typical super- 
novae in supcrgiant stars are not capable of cre- 
ating any relativistic ejecta (Matzner & McKcc 
1999), because there, shock acceleration ends at 
nonrelativistic velocities at a depth that matches 



the width of the shock. In this paper, we have 
only considered progenitors sufficiently compact 
and exploding with enough energy, that relativis- 
tic ejection does occur. 

Although we do not model asymmetrical explo- 
sions, our theory for shock and postshock dynam- 
ics can be applied to simulations of such an event. 
In §2.5 we present a formula (eq. [54]) for the rela- 
tivistic yield of each patch of the stellar surface, in 
terms of the intensity of the (nonrelativistic) shock 
as it might be observed in such a simulation. As 
multidimensional simulations cannot yet afford to 
treat relativistic dynamics with the extremely fine 
zoning that would be required to match the predic- 
tive power of this formula, it could be quite useful 
in future investigations. 

In §2.6 we have put our dynamical formulae to 
another use: the prediction of what will be ejected 
in the case of an explosion so weak that it only 
ejects a small mass from the outside of the star. 
Assuming ballistic motion after a brief phase of 
acceleration, the escape velocity marks the inner 
boundary of the ejecta. For this reason, the over- 
all kinetic energy (eq. [55]) and the kinetic energy 
above some given velocity (eqs. [56], [57a], and 
[57b]) are proportional the initial binding energy 
of the ejecta. Applied to neutron star phase tran- 
sitions in §3.3.2, equation (55) predicts within 50% 
the average kinetic energy per gram in the ejecta in 
Fryer & Woosley (1998)'s simulations. These con- 
siderations also provide some insight into the char- 
acteristic shock velocity at the base of the ejecta 
(§2.6.1). 

After calculating in §3.1.1 the energy distri- 
bution of relativistic ejecta from Woosley et al. 
(1999) 's model C06 for SN 1998bw, and demon- 
strating that this can explain both the observed 
gamma-ray burst (~ 35 s timcscale; §3.1.2) and 
the rapidly expanding radio shell (~ 12 day 
timescale; §3.1.3), we turned to the modeling of 
far more energetic "hypernova" explosions in §3.2. 
The association of GRB afterglows with the stel- 
lar population of their host galaxies (Bloom et al. 
2000b), and more tentatively with star forming 
regions (Kulkarni et al. 2000; Djorgovski, Bloom, 
& Kulkarni 2000; Galama & Wijers 2000), argues 
in favor of a connection with the deaths of massive 
stars and against scenarios involving neutron star 
mergers. Furthermore, the difficulty of avoiding 
baryon loading while trying to escape the stellar 



37 



progenitor, encountered even in highly asymmet- 
ric jet models (Aloy et al. 2000), is overcome in 
the hydrodynamic shock acceleration mechanism. 

In this model GRBs, like their afterglows, 
would represent radiation from an external shock 
that accelerates the ambient medium, with an ex- 
tra component (of comparable luminosity) from 
the reverse shock that decelerates the ejecta; both 
shocks must decelerate monotonically. Note that 
the higher efficiency of 7-ray production in mod- 
els including external shocks, compared to those 
based solely on internal shocks (Kumar 1999), 
gives external shock models less stringent total 
energy requirements. The spectral signature of 
such two-shock structures (Sari & Meszaros 2000) 
should be used to test this hypothesis; however, 
the reverse shock is buried under more material 
(by a factor of T) and may not be visible early 
on. While highly structured GRB profiles are dif- 
ficult to accommodate in this scenario (Sari & 
Piran 1997), most of the GRBs with afterglows 
have relatively smooth pulses. Small amounts 
of substructure may result from interaction with 
an inhomogeneous medium (Dermer & Mitman 
1999; Fenimore et al. 1999), or from postshock 
velocity perturbations produced by instabilities in 
the accelerating shock (Luo & Chevalier 1994) 
or the presence of steep density shelves in the 
progenitor's atmosphere. More detailed modeling 
of these processes is required. In some bursts, 
there is evidence that a smooth and soft compo- 
nent of the GRB flux is an early manifestation 
of the afterglow (Burenin et al. 1999; Giblin et 
al. 1999; Burenin 2000; Tkachenko et al. 2000). 
The existence of this component agrees with the 
external-shock model, although the fact that it 
does not account for the entire GRB does not. 
GRB 990123 had an optical flash that has been 
interpreted (e.g., Sari & Meszaros 2000) to arise 
from reverse-shock emission; however, this model 
invoked a brief rather than a continuous collision 
of ejecta with the swept-up shell. 

Because an external medium is necessary to lib- 
erate the energy, the Lorentz factor of the emitting 
region must not only exceed the usual lower limit 
from photon-photon opacity (eq. [69]), but also 
a limit due to the opacity of a circumstellar wind 
(eq. [71]) or the density of a uniform background 
(eq. [72] ) . The observed gamma-ray energies thus 
become, in a spherical model for the burst, lower 



limits for the amount of energy in ejecta whose 
Lorentz factor satisfies these constraints. As the 
bursts in question have quite smooth time profiles, 
and because our adoption of the external-shock 
hypothesis dictates that the age of a burst, rather 
than the timescale of its substructure, be used to 
constrain the Lorentz factor, we predict smaller 
minimum Lorentz factors than commonly quoted 
in the literature, by factors of ten or more. This 
eases the baryon-loading constraint by a similar 
amount. 

We find that a simple, and admittedly extreme, 
spherical hypernova model satisfies the energy and 
velocity constraints of many of the observed cos- 
mological GRBs. The most energetic bursts, such 
as GRB 990123, require asymmetric explosions. 
These need not be highly aspherical as the results 
of §2.5 show that a mild (factor of two) variation 
in shock intensity at a given column below the sur- 
face of the star can lead to an extreme (factor of 
100) variation in the observed burst intensity. If 
such variations are justified by simulations, this 
result may allow even the most energetic bursts 
to be explained by hydrodynamical shock accel- 
eration in aspherical hypernovae. Since the very 
luminous bursts do not differ from typical ones in 
any significant manner other than luminosity, it 
is likely that all cosmological bursts are aspher- 
ical, and as a result there is no need to invoke 
hypernova models as extreme as the one we have 
considered. 

Since we predict a kinetic energy that decreases 
with Lorentz factor roughly as 1/r^, there must 
be more energy in slower relativistic ejecta, and (if 
E- ln <C Mojc 2 ) much more in nonrelativistic ejecta 
as well. As the shocked region decelerates, an in- 
creasing fraction of this energy should appear in 
the afterglow; therefore, estimates of afterglow en- 
ergetics (e.g., Freedman & Waxman 1999) may 
potentially rule out such models in the future. 
However, having demonstrated in equation (73) 
that the energy of the shocked shell increases quite 
slowly, we find that Freedman & Waxman's con- 
straints cannot yet rule out hypernovae as sources 
of cosmological GRBs. 

Energetic events in compact objects have also 
been considered as engines for GRBs. White 
dwarfs (§3.3.1) explode as type la supernovae, but 
the relativistic ejecta from SNe la are too weak 
to be of much interest. Accretion- induced col- 
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lapse of white dwarfs can potentially concentrate a 
somewhat smaller energy (~ 10 50 erg) in a smaller 
(<~ 0.1 M ) mass of ejecta. Considering the sim- 
ulations of Fryer et al. (1999), we predict a much 
greater yield of relativistic ejecta than their es- 
timates; but not enough to power cosmological 
GRBs. 

Neutron star phase transitions, considered in 
§3.3.2, are well described by the weak explosion 
theory of §2.6. As equation (55) gives the overall 
energetics of the explosion to within 50%, we are 
justified in using equation (56) to predict the ener- 
getics of fast ejecta. Our estimates, validated with 
numerical simulation, are several orders of magni- 
tude higher than those of Fryer & Woosley (1998), 
who underestimated postshock acceleration. How- 
ever, we note that such events do not put enough 
energy at sufficient Lorentz factor to satisfy the 
observational constraints for distant bursts with 
afterglows. 

We thank Stan Woosley and Lee Lindblom for 
kindly providing models of SN 1998bw's progeni- 
tor star and a white dwarf, respectively. We also 
thank M. Aloy, A. Filippcnko, Z.-Y. Li and the ref- 
eree for detailed comments on the manuscript. We 
are grateful to A. Cumming, R. Fisher, P. Kumar, 
A. MacFadyen, R. Sari, E. Scannapieco, N. Sha- 
viv, and A. Youdin for helpful discussions. CDM 
appreciates the gracious hospitality of Sterl Phin- 
ney and Roger Blandford during his visits to Cal- 
tech, as well as financial support from NSERC. 
The research of JCT and CFM is supported by 
NSF grant AST 95-30480. 



REFERENCES 

Aloy, M. A., Miiller, E., Ibancz, J. M., Marti', J. M., & 
MacFadyen, A. 2000, ApJ, 531, L122 

Bailyn, C. D., Jain, R. K., Coppi, P. & Orosz, J. A. 1998, 
ApJ, 499, 367 

Baring, M. G. and Harding, A. K. 1997, ApJ, 491, 663 

Berezinsky, V. S., Blasi. P., &, Hnatyk, B. I. 1996, ApJ, 
469, 311 

Blandford, R. D., & McKee, C. F. 1976, Phys. of Fluids, 
19, 113 

Bloom, J. S. ct al. 1998a, ApJ, 506, L105 
Bloom, J. S. ct al. 1998b, ApJ, 508, L21 
Bloom, J. S. ct al. 1999, Nature, 401, 453 
Bloom, J. S. ct al. 2000, GCN 661 

Bloom, J. S., Kulkarni, S. R., & Djorgovski, S. G. 2000, 
astro-ph/0010176 

Briggs, M. S. et al. 1999, ApJ, 524, 82 

Burenin, R. A. et al. 1999, A&AS, 138, 443 

Burenin, R. A. 2000, Astronomy Letters, 26, 269 

Burrows, A. 1987, ApJ, 318, L57 

Castro, S. M. et al. 2000, GCN 605 

Ccntrclla, J. M. 1986, in Dynamical Spacetimes & Numer- 
ical Relativity, cd. Joan M. Centrella (Cambridge), 326 

Chandrasekhar, S. 1939, An Introduction to the Theory of 
Stellar Structure (New York: Dover) 

Chevalier, R. A. 1982, ApJ, 258, 790 

Chevalier, R. A. 1990, ApJ, 359, 463 

Chevalier, R. A., & Soker, N. 1989, ApJ, 341, 867 

Clayton, D. D. 1983, Principles of Stellar Evolution & Nu- 
cleosynthesis, Chicago 

Colgate, S. A. 1974, ApJ, 187, 333 

Costa, E. & Feroci, M. 1997, IAUC 6572 

Costa, E., Feroci, M, Piro, L., Soffitta, P., Amati, M., & 
Cinti, M. N. 1997, IAUC 6649 

Dadina, M. et al. 1999, IAUC 7160 

Dar, A., Kozlovsky, B., Nussinov, S., & Ramaty, R. 1992, 
ApJ, 388, 164 

Dcrmcr, C. D., & Mitman, K. E., 1999, ApJ, 513, L5 

Dingus, B. L. et al. 1997, Gamma-Ray Bursts: J^th 
Huntsville Symposium, eds. Charles A. Mccgan, Robert 
D. Preece & Thomas M. Koshut. Woodbury, New York 
: AIP 428 



39 



Djorgovski, S. G., Kulkarni, S. R., & Bloom, J. S. 1999, 
GCN 289 

Djorgovski, S. G. ct al. 1999, GCN 481 

Djorgovski, S. G., Bloom, J. S., & Kulkarni, S. R. 2000, 
ApJ, submitted (astro-ph/0008029) 

Esin, A. A. & Blanford, R. D. 2000, ApJ, 534, L151 

Fcnimore, E. E., Cooper, C, Ramirez-Ruiz, E., Sumner, 
M. C, Yoshida, A., & Namiki, M. 1999, ApJ, 512, 683 

Fcroci, M. ct al. 1999, IAUC 7095 

Freedman, D. & Waxman, E. 2000, ApJ, submitted, astro- 
ph/9912214 

Fruchter, A. S. et al. 1999, ApJ, 516, 683 

Fryer, C, & Woosley, S. E. ApJ, 501, 780 

Fryer, C, Benz, W., Herant, M., & Colgate, S. A. ApJ, 
516, 892 

Fynbo, J. P. U. et al. 2000, GCN 807 

Galama, T. J. et al. 1998, Nature, 395, 670 

Galama, T. J. ct al. 2000, ApJ, 536, 185 

Galama, T. J., & Wijcrs, R. A. M. J. 2000, ApJ, submitted 

Germany, L. M. et al., 2000, ApJ, 533, 320 

Giblin, T. W., van Paradijs, J., Kouvcliotou, C, Con- 
naughton, V., Wijcrs, R. A. M. J., Briggs, M. S., Preece, 
R. D. & Fishman, G. J. 1999, ApJ, 524, L47 

Glcndcnning, N. K. 1997, Compact Stars (Springer, New 
York) 

Gnatyk, B. I. 1985, Sov. Astron. Lett. 11(5), 331 

Goodman, J., Dar, A., & Nussinov, S. 1987, ApJ, 314, L51 

Hamann, W.-R., Koesterke, L., & Wessolowski, U. 1995, 
A&A, 299, 151 

Harpaz, A. 1994, Stellar Evolution, Peters. 

Harrison, F. A. et al. 1999, ApJ, 523, L121 

Heger, A., & Langcr, N. 1996, A&A, 315, 421 

Heise, J. & in 't Zand, J. 1997, IAUC 6787 

Hoflich, P., Wheeler, J. C, & Wang, L. 1999, ApJ, 521, 
179 

Hoflich, P., Wheeler, J. C, Hincs, D. C. & Trammell, S. R. 
1996, ApJ, 459, 307 

Hogg, D. W. & Fruchter, A. S. 1999, ApJ, 520, 54 

Hurley, K. et al. 1994, Nature, 372, 652 

Hurley, K. & Cline, T. 1999, GCN 450 



Hurley, K., Cline, T., & Mazets, E. 2000, GCN 642 

Hurley, K., Mazets, E., Golcnctskii, S., & Cline, T. 2000, 
GCN 801 

Iwamoto, K. ct al. 1998, Nature, 395, 672 

Johnson, M. H., & McKee, C. F. 1971, PRD, 3, 858 

Kazhdan, Y. M., & Murzina, M. 1992, ApJ, 400, 192 

Khokhlov, A. M., Hoflich, P. A., Oran, E. S., Wheeler, J. 
C, Wang, L. & Chtchelkanova, A. Y. 1999, ApJ, 524, 
L107 

Koesterke, L., & Hamann, W. -R. 1995, a, 299, 503 

Krolik, J. H. and Pier, E. A. 1991, ApJ, 373, 277 

Kulkarni, S. R. et al. 1998, Nature, 395, 663 

Kulkarni, S. R. et al. 1999, Nature, 398, 389 

Kulkarni, S. R. ct al. 2000, in Gamma-Ray Bursts: 5th 
Huntsville Symposium, astro-ph/0002168 

Kumar, P. 1999, ApJ, 523, L113 

Li, Z. Y., & Chevalier, R. A. 1999, ApJ, 526, 716 

Lindblom, 1999, Phys. Rev. D, 60, 064007 

Lithwick, Y., & Sari, R. 2001, ApJ, submitted, astro- 
ph/0011508 

Litvinova, I. Y. & Nadfizhin, D. K. 1990, Soviet Astron. 
Lett., 16 29 

Luo, D. & Chevalier, R. A. 1994, ApJ, 435, 815 

MacDonald, D. A., Thorne, K. S., Price, R. H., & Zhang, 
X.-H. 1986, in Black Holes, the Membrane Paradigm, 
cd. K. S. Thorne, R. H. Price, & D. A. MacDonald 
(New Haven: Yale Univ. Press) 

MacFadycn, A., & Woosley, S. E. 1999, ApJ, 524, 262 

MacFadyen, A., Woosley, S. E., & Heger, A. 2000, ApJ, 
submitted, astro-ph/9910034 

Marcaide, J. M. and 15 colleagues 1997, ApJ, 486, L31 

Marti, J. M., & Miiller, E. 1994, J. Fluid Mech., 258, 317 

Matzner, C. D., & McKcc, C. F. 1999, ApJ, 510, 379 

Meszaros, P. & Rees, M. J. 1993, ApJ, 405, 278 

Meszaros, P. 1999, Nature, 398, 368 

Nakamura, T., Mazzali, P. A., Nomoto, K., & Iwamoto, K. 
2000, ApJ, submitted 

Nadyozhin, D. K. 1985, Ap&SS, 112, 225 

Norris, J. P., Bonncll, J. T. & Watanabe, K. 1999, ApJ, 
518, 901 



40 



Norris, J. P., Marani, G. F. & Bonnell, J. T. 2000, ApJ, 
534, 248 

Ostriker, J. P., & McKee, C. F. 1988, Rev. Mod. Phys., 60, 
1 

Paczynski, B. 1998, Gamma-Ray Bursts: 4th Huntsville 
Symposium, cds. C.A. Meegan, R.D. Preece, & T.M. 
Koshut, Woodbury: New York 

Panaitescu, A., Meszaros, P. & Rees, M. J. 1998, ApJ, 503, 
314 

Patat, F., & Picmonte, A. 1998 IAU Circ. 6918 

Pendleton, et al. 1997, Gamma-Ray Bursts: 4th Huntsville 
Symposium, cds. C.A. Meegan, R.D. Preece, & T.M. 
Koshut, Woodbury: New York 

Pian, E. & 27 colleagues 2000, ApJ, 536, 778 

Piran, T. 1999, Physics Reports, 314, 575 

Rees, M. J. & Meszaros, P. 1998, ApJ, 496, LI 

Reichart, D. E. 1999, ApJ, 521, 111 

Sakurai, A. 1960, Comm. Pure Appl. Math., 13, 353 

Sari, R., & Piran, T. 1997, MNRAS, 287, 110 

Sari, R. & Meszaros, P. 2000, ApJ, 535, L33 

Schneid, E. J. et al. 1995, ApJ, 453, 95 

Sedov, L. I. 1946, Prikl. Mat. Mekh. 10, 241, No. 2 

Sedov, L. I. 1959, Similarity & Dimensional Methods in 
Mechanics (Academic, New York) 

Shapiro, P. R. 1979, ApJ, 233, 831 

Shaviv, N. J. & Dar, A. 1995, ApJ, 447, 863 

Smith, D. A., Hurley, K., & Cline, T. 2000, GCN 568 

Soffitta, P., et al. 1998, Mf7C6884 

Sollerman, J. et al. 2000, ApJ, 537, L127 

Springmann, U. 1994, A&A, 289, 505 

Taylor, G. I. 1950, Proc. R. Soc. London, 201, 159 

Thorsett, S. E., & Hogg, D. W. 1999, GCN 197 

Tinney, C, Stathakis, R., Cannon, R., & Galama, T.J. 
1998, IAUC 6896 

Tkachcnko, A. Y., Tcrckhov, O. V., Sunyaev, R. A., Bu- 
renin, R. A., Barat, C, Dezalay, J. -. & Vedrcnnc, G. 
2000, A&A, 358, L41 

Turatto, M. ct al. 2000, ApJ, 534, L57 

Wang, L. & Wheeler, J. C. 1998, ApJ, 504, L87 

Wang, L., Howell, D. A., Hocflich, P., & Wheeler, J. C. 
2000, ApJ, submitted, astro-ph/9912033 



Waxman, E., & Loeb, A. 1999, ApJ, 515, 721 
Waxman, E., & Draine, B. T. 2000, ApJ, 537, 796 
Weaver, T. A. 1976, ApJS, 32, 233 

Weiler, K. W., Panagia, N., Sramek, R. A., Van Dyk, 
S. D., Montes, M. J., & Lacey, C. K. 2000, in The 
Largest Explosions Since the Big Bang: Supernovae 
and Gamma Ray Bursts, 1999 STScI May Symposium, 
astro-ph/0002501 

Wheeler, J. C, Yi, I., Hocflich, P., & Wang, L. 2000, ApJ, 
537, 810 

Wheeler, J. C, Hoeflich, P., Wang, L., & Yi, I. 2000, in 
Gamma-Ray Bursts: 5th Huntsville Symposium, astro- 
ph/9912080 

Wieringa, M. H., Kulkarni, S. R. & Frail, D. A. 1999, 
A&AS, 138, 467 

Woosley, S. E. 1993, ApJ, 405, 273 

Wooslcy, S. E., & Weaver, T. A. 1986, ARA&A, 24, 205 

Woosley, S. E., & Baron, E. 1992, ApJ, 391, 228 

Woosley, S. E., Eastman, R. G., & Schmidt, B. P. 1999, 
ApJ, 516, 788 



This 2-column preprint was prepared with the AAS IATjrjX 
macros v5.0. 



41 



